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ABSTRACT 


This thesis investigates various techniques for fore- 
easting the volume of containerizable cargo that flows into 
the container stuffing station at the Military Ocean 
Demminal, Bay Area, Oakland, California. Cargo input data 
Bemanalyzed in terms of weekly cargo volume inputs for a 
selected number of major ports of debarkation. The time- 
series data for these POmseomrsSatarst tested for serial 
Semrelation. Based on the affirmative results of the 
serial correlation test, the following forecasting methods 
Pre investigated: the moving average, the exponentially 
weighted average, the exponentially weighted average with 
trend adjustment and the exponentially weighted average with 
an adaptive response rate. By means of statistical testing 


procedures, the "best" forecasting method is determined. 
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I. INTRODUCTION 


imimeuaiacect the dispersal of United States forces 
throughout the world, the shipment of the vast quantities 
of needed supplies and material to overseas stations is a 
predominant logistics Problcimnonme the Department of Defenge. 
Since the late 1960's the need for efficient and economical 
manmdling of this cargo has led the DoD to use containerized 
‘shipments Simeone teral=Sea-POIng Carriers as the primary 
fotane transportation source. The retention of break bulk 
type shipping operations has been necessary due to the size 
imetations of containers. Also, the fulfillment of urgent 
memlmrements where time 1S critical have been met using air 
Mmeamsport. However, ocean containerization offers the most 
@est-etfective method of handling routine shipments between 
Pemts equipped to handle this type of traffic, as are most 
metric Ocean ports. 

Miewabality to fill or stuff a container at the con- : 
Seenon's warehouse and unpack it at the consignee's ware- 
house, with no intermediate handling of the cargo, is a 
tremendous asset. However, many shipments are not large 
enough to economically warrant their own vans; therefore, 
stuffing with other shipments bound for the same destination 
may be accomplished at a central location. One such stuff- 
PMipetdctlityess the container freight division (CFD) of the 
Military Ocean Terminal Bay Area (MOTBA) located in Oakland, 


California. This facility stuffs vans with containerizable 
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Cargo enroute to ports, mainly in the Pacific and the Far 
Paeeeeand arranges transportation of the containers to the 
Pomes Via commercial carrier. 

All commercial carriers require that bookings for vans 
and space aboard vessels be made in advance of the ship's 
loading date. As of the fall of 1974, this lead time was 
at least three weeks. To effectively make these bookings 
memane advanced date, a prediction of the cargo on hand on 
the date of arrival of the vans is required. Normally, 
vans arrive at the stuffing station four to six days prior 
to the loading date. Thus, the request for space must take 
meer account {1) the present volume of cargo in the ware- 
house, (2) the forecasted arrivals of cargo during the 
famleanpecween che booking and the van arrivals for stuff- 
ing and (3) any previous bookings that will deplete the 
eareo volume destined for the port during the forecasting 
period. 

Piemzeneral objective of this thesis is to examine 
practical forecasting methods which can be utilized to pre- 
dict cargo arrivals at the container stuffing station. To 
predict arrivals there are two general approaches that 
Semlanbe used. The first approach is based on notification 
of shipment departures by consignors, which would provide 
Memeaiced Notice of arrivals at the stuffing station. The 
Other approach is based on historical data, from which pro- 
Meetions for future arrivals may be drawn. 

An advance notification system does not presently exist 


for the majority of cargo arriving at the station; however, 
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it could be produced by requiring consignors to provide the 
CFD with three-week advance notification of cargo shipments, 
perhaps by post card. Use of this approach would require 
analysis of the delay between the arrival of notifications 
and shipments from each consignor. Also, the statistical 
@amacteristics and stability of the transit times for ship- 
Ments and notifications must be determined. It is doubtful 
that most consignors have enough advance knowledge of a 
shipment's departure to provide a timely and accurate noti- 
fication, especially if the departure is two or three weeks 
in the future. The effect of the possible lack of consignor 
cooperation must also be assessed. Taking into account the 
mesources that would be required to effectively implement 

an advance notification approach, it was decided that fore- 
casting methods which rely only on historical data for 
future projections would be more appropriate for a thesis 
topic. 

The historical data used here consisted of the cargo 
arrival data for the CFD for fiscal year 1974; it repre- 
sented the most recent peacetime data available. The data 
provided the volume, weight and date of arrival of each 
shipment. Of the three alternative units for forecasting 
(number of shipments, volume and weight), the volume of the 
mensO Was selected because it is currently used by the mari- 
time community as the measure for booking containers on 
messels. 

Two alternative aggregation levels could have been used 


for forecasting. One was based on the port of debarkation 
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(POD), and the other based on consignee destination of the 
@ameo. Lhere are more than 1,000 active consignees served 
by the CFD, which would require a very large number of in- 
dividual forecasts if the consignee base were used. PODs 
are hie currently used aggregation levels for forecasts by 
CFD personnel because van space is, by convention, booked 
aboard ships destined for PODs, not consignees. Therefore, 
PODs were selected as the aggregation level. Analysis of 
the data revealed 16 major PODs which received over 90% of 
the shipments passing through the CFD. 

The effectiveness criteria used for the selection of an 
optimal method are: 1) computational feasibility, and 2) 
the ability of the method to maximize forecasting accuracy 
by minimizing forecasting error. Computational feasibility 
1s limited here to those methods which do not require the 
use of a computer to routinely generate forecasts. Mean 
absolute deviation and mean deviation squared are used as 
Measures ot forecast errors. 

Forecast errors were computed over time horizons of 
one, two and three weeks, corresponding to possible lead 
Mees in the booking process. In other words, for a two- 
week time horizon a single forecast was used for two suc- 
cessive weeks. The deviation for the forecast was computed 
by taking the average difference between the forecast and 
the actual input volume for each of the weeks covered by 
the forecast. 

Generally, the selection among various methods of fore- 


Sa-iieeaepends Upon the statistical characteristics of the 
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meee these Characteristics are differentiated by the degree 
Mamsiatistical dependence between successive observations. 

If there is no statistical dependence between successive 
observations, then ordinary trending and indexing techniques 
meaepe applied. On the other hand, the presence of statis- 
tical dependence among successive observations, (called 
serial correlation), implies a need for methods which place 
more emphasis on recent past observations when predicting 
future volumes. 

Tests for serial correlation of the data for the major 
Pelsewere made using the Durbin-Watson d-statistic. The 
memereoated results of these tests indicated that the pre- 
mamtecnot scrial correlation was statistically significant. 

Figure 1 is a graph of the weekly input volumes of POD 
Meeetor a, 20-week period commencing with the nineteenth 
week of FY '74. This was selected as a typical representa- 
meomeot the data for the major PODs. This particular POD 
exhibited strong positive serial correlation. 

Wieieune Existence of serial correlation in the data 
confirmed, the moving average technique and three exponen- 
tially weighted average techniques were examined for pre- 
diction accuracy. These methods were chosen for their 
Beobative ease Of Computation and ability to adapt to changes 
in volume inputs over time. 

Simple moving averages [Ref. 5] of from one to ten week 
durations were the first methods explored. In each method, 
the most recent weeks for the designated duration are aver- 


aged and the result applied to the forecast through the 
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Figure 1: Weekly volume inputs for POD UL/ 
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designated horizon. The three-week moving average was the 
technique in actual use at the CFD. The resulting data in- 
dicated that for all time horizons, moving averages greater 
than three weeks generally provided lower deviations and 
thus, better forecasts than the three-week moving average. 
Note that the longer the duration of the average, the less 
meeponsive the predictions are to recent changes. The data 
indicated that over-sensitivity was as undesirable as no 
sensitivity in most cases. 

iiemnext method examined was the simple exponentially 
weighted average, [Ref. 5], or exponential smoothing method. 
This method used a weighting or smoothing constant, alpha, 
@ Value usually between zero and one, to regulate the em- 
phasis on the most recent observation. Large vaiues of aipha 
place increasing emphasis on the most recent observations 
and less on preceding observations. A forecasted average 
for an observation period was computed as the sum of the 
observations in that period times alpha, plus the forecasted 
average for the previous period times one minus alpha. In 
this way, the new observation was added to the accumulated 
weighted average, a continually changing number. This fore- 
cast average was then used as the forecast and applied over 
the next weeks of the designated horizon and the deviation 
Beatistics accumulated. 

To discover the levels of alpha which minimized error 
measurement, the levels of alphas were varied by 0.01 from 
auee tO 1.0, and applied to the data for each POD and hori- 


zon. The alpha levels which produced minimal deviation 
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meristics ranged from -0.0Z2 to +0.77 depending on the POD 
and time horizon specified. The optimal single alpha level 
for all PODs was about 0.35 and varied slightly depending 
on the time horizon. 

The next forecasting method explored was the exponen- 
tially weighted average with adjustment for trends [Ref. 5]. 
Repco procedure was examined to determine if trend adjust- 
Menmtes would reduce forecast error. This procedure started 
with the same basic formula as the simple exponentially 
weighted average, however a trend factor is also computed. 
This trend is accumulated and smoothed in the same manner 
as was applied to the forecast averages, resulting in an 
average trend adjustment. The forecast volume to be used 
mimeme yredictions was created by combining the forecasted 
average with the average trend adjustment. 

To obtain minimal errors, alpha levels were again varied 
Seve ol trom -0.10 to 1.00 and the method applied to the 
data for all PODs. The results revealed minimized devia- 
tions for alpha levels between -0.01 and 0.38, again, de- 
Mending on the POD and time horizon specified. When 
accounting for trends, the optimal single alpha level for 
metros was 0.16, varying slightly for the different time 
Morizons. 

The final forecasting procedure used was the exponen- 
tially weighted average using an adaptive response rate 
technique to modify the alpha level used from period to 
period to account for fluctuations in the observations 


[Ref. 18]. The absolute and actual errors, or deviations 


19 





and actual errors between the observed and forecasted 

volumes were separately accumulated and smoothed in the 

same fashion as the forecast average using a smoothing con- 
stant, gamma. A tracking signal was then computed by di- 
viding the smoothed absolute error by the smoothed error. 
Deviation statistics using the simple exponentially weighted 
average were computed with the alpha level for each forecast 
set at the absolute value of the tracking signal. 

To evaluate the adaptive response rate method, it was 
em@erred to the data for each of the PODs while varying the 
time horizon and the level of the smoothing constant, gamma, 
by 0.01 from -1.00 to 1.00. Here the use of negative con- 
stants indicates that the new smoothed error term is created 
Peeenlarging the ola smoothed error and adding or subtract- 
ing a portion of the new error to it. The minimum deviation 
Statistics for the individual POD and various time horizons 
were realized with gamma levels ranging fon porate Ole 
The resultant optimal single gamma level for all POD was 
about -0.05 depending on the horizon used. 

Mmecouwimiary, the experimentation with the above forecast- 
ing methods had created deviation statistics for the follow- 
ing forecasting procedures: 

1) Three-week moving averages (maintenance of status 
quo) ; 

2) Moving averages with a single, aggregately optimal 
mirat ion : 

3) Moving averages with the optimal duration for each 


Pop: 
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4) Simple exponentially weighted averages using a 
single, aggregately optimal alpha level; 

5S) Simple exponentially weighted averages using in- 
dividually optimal alpha levels for each POD; 

6) Exponentially we1ghted averages with trend adjust- 
ment using a single, aggregately optimal alpha level; 

io Exponentially medeaced averages with trend adjust- 
ment using individually optimal alpha levels for each POD; 

8) Exponentially weighted averages with adaptive 
response rates using a single, aggregately optimal ee 
level; 

9) Exponentially weighted averages with adaptive re- 
Sponse rates using individually optimal gamma levels for 
each POD. 
eeme these deviation statistics, the selection of the "best" 
forecasting method was made subject to the following 
considerations: 

1) For operational consistency, it was assumed that 
only one of the methods would be chosen for use with all 
PODs. No attempt was made to evaluate the problem by 
choosing the optimal forecasting method for each POD. 

2) The deviation statistics for the one week time 
horizon were used in the evaluation as the relative effec- 
mveness of the various methods for each POD. 

3) The mean absolute deviation statistic was selected 
feetede Principal goodness criterion since no extra penal- 


ties for large individual forecast errors were intended. 
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imitenewend, statistical significance tests are used to 
reveal whether the mean absolute deviations using one fore- 
casting method were significantly less than the deviations 
of another method. The tests were conducted at a signifi- 
cance level of 0.05. The results indicated that the simple 
exponentially weighted average and adaptive response rate 
methods were not significantly different from each other. 
However, both methods were better than any of the others. 

Based On the above results the use of the simple expo- 
nentially weighted average was chosen. The selection was 
moumaed On the ease of computation of the forecast, as the 
adaptive response rate method required many more steps for 
each calculation and offered no better results. Addition- 
ally, the adaptive response rate method required a broad 
range of gamma levels and was characterized by multi-modal 
relationships between the gamma level and MAD for most cases 
of PODs and time horizons. This indicated great sensitivity 
of accuracy to the gamma levels for particular data sets. 

When compared to the three-week moving average, the use 
of this simple exponentially weighted average technique 
provided 4.9%, 8.0% and 10.0% reductions in mean absolute 
deviations for the one-, two- and three-week time horizons 
respectively. In addition to increasing the forecast ac- 
curacy the exponentially weighted average is also as easy 
computationally as the three-week moving average. Opera- 
Buona Denerits of the increased forecast accuracy are a 
reduction in the average age of cargo and volume of cargo in 


the warehouse, and increased cube utilization. 
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II. BACKGROUND 


A. THE BEGINNING OF CONTAINERIZATION 

While containers had been used for transporting cargo 
aboard vessels as early as the 1820's, the modern era of 
containerization began in 1966 when Malcolm McLean of Sea- 
Land Service, Inc., began loading containers aboard ships 
moma £Or Europe. The economics of contalinerization soon 
became apparent. By the end of 1968, one-fourth to one- 
third of shipments on major U.S. routes were containerized. 
American-flag operators as well as their foreign-flag com- 
petitors began to invest heavily in the containerized sys- 
meme this system includes not only containers, but also 
Specialized ships, container terminals and road equipment. 

Within the industry the container came to be called the 
greatest thing in packaging since the paper bag. Most con- 
tainers are simple boxes of plywood or aluminum and are 
Memm@eactured in standard sizes: approximately 8 feet wide, 
8 feet deep and various lengths between 20 to 40 feet long. 
The container volume capacity ranges from 1,100 cubic feet 
mom ,t00 Cubic feet, and weight capacity ranges from 38,000 
pounds to 46,000 pounds depending on the size of the van. 

Containerization required the shipping industry to 
adopt an entirely new line of CrbiciiIomana provided the in- 
dustry with many phases of economic improvement: more ef- 
Peeteney sin careo and terminal handling, a reduction in 


merchandise damage due to weather and handling, a reduction 
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in terminal facility requirements, a savings in reduced 
special packaging for export, a reduction in pilferage and 
other losses, the ability of manufacturers to operate with 
meameed inventory, and the benefit of computerized inven- 
Memy=-CcOntrol systems. 

Containerization has cut shiploading and unloading time 
from as much as eight days to as little as twelve hours. 
The varied forms of containerization are enabling cargo 
ships to spend as many as nine out of every ten days' pro- 
Mmetmivity at sea instead of half their time in port -- 
Causing some shipping lines to have alternating crews so 
that the men can spend some time with their families. A 
typical containership is loaded with some 700 vans, seven 
deep and six across, and is on its way at 22 knots. 

The new era has resulted in transmodal systems to move 
Sareo Via land, sea and air with computerized efficiency. 
In some cases, overall transport costs are down 30-40% as 
compared to conventional break bulk cargo handling. 

Ped=elband ocrvice, Inc., has become the biggest container 
operator in the world and is one of the two largest U.S. 
steamship lines. Sea-Land has five 33-knot containerships 
With capacity for 1,082 units. With other carriers taking 
advantage of container economics, the United States flag 
Simrters Nave become the world's largest containership 
fleet. 

The port of Oakland, California is now the leading West 
Coast cargo shipping terminal -- mainly due to containeriza- 


tion. The port handles 45% of all container shipping on 
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the West Coast and is the second largest container port in 
the world -- New York being the largest. 

According to James L. Reynolds, President of the 
pi@emican Institute of Merchant Shipping (AIMS), by 1976 
Sameo moving to and from U.S. ports is expected to rise to 
270 million measurement tons (a unit of volume representing 
fmemprc feet) plus 7.5 million in military cargo and an- 
other 20.5 million in government impelled cargo such as 
foreign aid shipments. 

Containerization has been recognized as a system that 
will best serve the needs of the puUolEe=andethe= transport a- 
tion industry. The container revolution stems from the im- 
mense transportation capability gained from moving goods 


swiftly, safely and reliably. 


B. CONTAINERIZATION AND THE DEPARTMENT OF DEFENSE 

mumucecntnyecarsecne United States Department of Defense 
has used commercial containerization as the primary means 
of transporting its general ocean cargo. Substantial cost 
aviiies trom the use of containerships in lieu of using the 
traditional break bulk cargo ships have been recognized. 
Memeencre Still remain problems for transportation offi- 
Sials within DoD to reckon with. 

PpeuGeezeitiustrates the portion of the transportation 
SS eccm Within the Department of Defense to which this thesis 
addresses itself. 

Orders are placed by DoD overseas commands known as 


consignees. Various vendors and military depots known as 
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consignors, fill these purchase orders which the consignees 
Meea to execute their missions. It now becomes the job of 
the DoD transportation managers to deliver these items to 
the consignee. 

Cargo leaving the consignor for the consignee may be 
Peamsported across the ocean in either container or break 
bulk ships. Conventional ships transport break bulk cargo 
to the requesting overseas command. If the goods are 
Ordered in large enough quantities and are containerizable, 
the consignor may place the cargo in containers at his ware- 
house. This is known as source stuffing. The sealed con- 
tainer is then shipped to the POE where it has been booked 
aboard a containership for overseas transport to the con- 
Signee. This type of operation provides the most rapid and 
economical form of ocean shipping. 

Another situation occurs when container stuffing sta- 
memos, at the POE or inland, receive quantities of contain- 
erizable cargo for consolidated stuffing into containers. 
At the stuffing station, cargo for particular consignees is 
stuffed only after meeting volume and weight restrictions. 
A commercial shipping company then receives the containers 
meomethe stuffing station for lift aboard containerships 
mcm overseas transport. 

Pemthe Port of Debrakation (2OUiecinglcweonsitpenee Ccon- 
tainers are sent directly to the respective consignees. 
Containers with cargo for more than one consignee are trans- 
merted to designated break bulk stations (BBS), which nee 


itself by a consignee, where they are brokem down for further 
shipment to individual consignees. 
Le 





C. WESTERN AREA, MILITARY TRAFFIC MANAGEMENT COMMAND 

Western Area, Military Traffic Management Command (WAMTMC) 
Hema jOintly-staffed field organization under Headquarters, 
Military Traffic Management (MTMC), Washington, D. C. WAMTMC 
is responsible for transporting export cargo of the Depart- 
ment of Defense and other government agencies to overseas 
commands. It has been located at the Oakland Army Base, 
Oakland, California, since being established in February 
(1965. MTMC is a major field command of the Department of 
the Army. 

WAMTMC is responsible for transportation management of 
domestic and export shipments in 14 western states: Idaho, 
Arizona, California, Montana, Nevada, Oregon, Utah, New 
Mexico, Washington, Colorado, Nebraska, South Dakota, North 
Dakota and Wyoming. One of the three Military Ocean Termi- 
Nals Operated by WAMTMC on the West Coast will receive ex- 
bore cargo from these 14 states which is destined for 
shipment to overseas installations. The three terminals 
m@eemated by WAMIMC are: Military Terminal Unit, Pacific 
Northwest, Seattle, Washington (PNW); Military Ocean Termi- 
Males Bay Area, Oakland, California (MOTBA); Southern 
fertormia Outport, Long Beach, California (SCO). During 
peacetime, WAMTMC also directs the flow of all air cargo 
in the Military Airlift Command's transportation system 


@mrouch the Military Airlift Clearance Authority (MACA). 


Pee MILITARY OCEAN TERMINAL, BAY AREA 
Military Ocean Terminal, Bay Area is the consolidation 


of the Army and Navy terminal facilities in the San Francisco 
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Pewee a. it represents the largest operating unit of 
WAMTMC with terminal facilities at the Oakland Army Base 
Sa the Alameda Reefer Facility. Seven deep water berths 
are under its control as well as an 84-acre tidewater con- 
tainer stuffing area located at the Oakland Army Base. 
Also under MOTBA's control are the ports of Stockton, 


Sacramento, and Eureka. 


E. TIDEWATER CONTAINER STUFFING STATION 

The Container Stuffing Station (CSS) operates under the 
Gontainer Freight Division (CFD) at MOTBA. Containerizable 
cargo is loaded into the shipping containers at both instal- 
Hatzons. These containers are then transferred to the com- 
mercial shipping companies for shipment to overseas commands. 

iiewGos 1S Operated under civilian contract. Container- 
izable cargo is placed in warehouses until containers become 
available for stuffing. In most cases commercial containers 
are used for stuffing, although military vans are sometimes 
Meeds the contractor has certain restrictions on container 
Seurting, some examples are cargo mix, container utiliza- 
muon, and consignee mixing eemcacions . Mimadd ton Velume. 
weight and size of containers vary among the various com- 
Gememal Carriers. There are 16 types of containers presently 
in use; and, in general, each company can accommodate only 
those containers which have been specifically designed for 


its vessels. 
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F. MILITARY SEALIFT COMMAND, PACIFIC 

The Military Sealift Command »(MSC) 1s an agency that 
coordinates the various ocean transportation activities of 
the Department of Defense. Within the continental United 
States (CONUS) the relationship between the Military Sealift 
Command, Pacific (MSCPAC) and WAMTMC is very close. In 
Order to get its cargo aboard break bulk or container ves- 
sels, WAMTMC notifies MSCPAC of the amount of space required, 
the number of containers needed by the shipper and the size 
of the containers. MSCPAC books the space aboard a vessel 
and notifies WAMTMC of the vessel sailing date, lift date 
for the containers and the number of containers booked. 
This procedure is carried out whether the shipper is a ven- 


dor or WAMTMC is the shipper. 


G. MTMC MANAGEMENT INFORMATION SYSTEMS 

MTMC employs two management information systems in order 
memcontrol the flow of cargo within CONUS that will be 
eventually shipped to overseas commands. Both systems 
Moemerze a Burroughs 5500 computer. The two systems are the 
Meeianized Export Traffic System (METS) and the Surface Ex- 
port Cargo System (SURS). 

Heleemontcors the Bol’ cargo flow from its origination 
at the consignor's warehouse until its arrival at a state- 
Side Port of Embarkation, i.e., one of WAMTMC's three ocean 
terminals. The consignor, which may be a commercial vendor 
Tectibiicary depot, thus provides the basic input into the 


System. This input includes shipment data, commodity 
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meer mation, transportation data, destination, etc. Im- 
Beemerly formated documentation is frequently received from 
the commercial vendors requiring the correction of these 
errors by manual methods. 

Once cargo is within one of the three west coast ocean 
terminals, it falls under SURS, or SURS/CARDPAC. A CARDPAC 
is a set of eight IBM cards containing various shipping 
data. These cards have been punched by terminal personnel 
with data obtained from METS and any additional information 
that may be needed from shipping forms. 

PMtcmGnikNPAG; the Expected Receipt List (ERL), a list 
Meche Cargo expected in a particular shipment; and the 
Meeercable source document, ¢.g., a Government Bill of 
Meaane (GBL), Commercial Bill of Lading (CBL), Dray Tag, 
or Transportation Control and Movement Document (TCMD) are 
Boawarded to the terminal's cargo receiving area. As the 
Cargo is unloaded a check is made against information con- 
tained on the CARDPAC and corrections are made if necessary. 
The CARDPAC is then attached to the shipment and will be 
Meatized to monitor its flow through the terminal facilities. 

Mimotimmary, by using the METS and SURS management infor- 
mation systems, WAMTMC is able to monitor the flow of cargo 
into and within its terminal facilities. Management reports 
generated by the system and statistical analyses based on 
the information provided by the system are only as good as 


mie data going into it. 
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fee CARGO CLASSIFICATION 

There are two types of cargo passing through the Con- 
Geaaner Stuffing Station: (1) release unit (RU); and (2) 
less than release unit (LRU). 

Shipments in excess of 10,000 pounds and special cate- 
mermyecargo, such as classified cargo, are classified as 
Bemease Unit material. Under DoD Regulation 4500.32-R, 
Military Transportation and Movement Procedures (MILSTAMP), 
mmease Unit material requires positive export traffic 
mebease. 

Consignors having break bulk or containerizable cargo 
Maat) tall into RU classification must request and receive 
Srearnace to move that material before it can be shipped. 
WAMTMC, Export Control Division, serves as the clearance 
menor ity for RU material for the fourteen Western states. 

This thesis is not concerned with break bulk or source 
Seurted cargo. Therefore, this paper will not go into a 
detailed description of the booking procedure except to 
femeron that 1t iS accomplished through the coordinated ef- 
forts of MSCPAC, WAMTMC, the shipper and the commercial 
Sarrier. 

meesnould be mentioned that 1f RU classified cargo is 
itwemot SULfiCieént quantity to be source stuffed, it is 
emipped to the CFD for stuffing. 

Meityepercent Of the shipments passing through MOTBA 
are of the LRU classification. LRU material does not re- 
quire positive release into the system. Therefore, there 


1S an uncontrolled flow of cargo coming into the Container 
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Freight Division. This results in a very difficult booking 
problem for the Container Freight Division, MOTBA, WAMTMC 
ama@mmMocPAC, it is the forecast of this uncontrolled cargo 
flow to which this thesis addresses itself. 
pewircmomyepresciits the Cango Cycle to be considered. 
Meemperecedures at each step are as follows [Ref. 15]: 

1) <A DoD overseas command makes a request for goods 
Dy requisition. The specific requestor is the user of these 
goods and the ultimate consignee. The requisition is to be 
filled by the shipper (consignor) who may be a commercial 
vendor or a military depot. | 

Pee Material that 1s classified as less than release 
unit is shipped by the consignor directly to the Container 
Peemeht Division at MOTBA. 

3) Under Military Traffic Management Regulation, con- 
tainer requirements are offered to WAMTMC by MOTBA. 

meeeerequest for booking is then submitted to MSCPAC. 

eye MSCPAC then books space aboard ocean carrier. 

6) Once booking aboard the vessel is made, the related 
information is forwarded to WAMTMC. 

7) The Container Freight Division at MOTBA receives 
the booking data from WAMTMC. 

8) The Container Freight Division communicates directly 
with the ocean carrier to arrange for spotting and pick-up 
Of containers. 

9) Once containers have been stuffed at the CSS, they 
macecelivered directly to the ocean carrier's container yard. 


Promethis point shipment from POE to POD is made. 
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10) Delivery of material is made to the consignee by 
either the ocean carrier or by arrangements made by the 


overseas command. 


I. CURRENT BOOKING PROCEDURES 

As previously stated, a request for booking is made 
twenty-one days in advance and is based on a forecast of 
the volume of cargo in measurement tons that is expected to 
be On-hand at the time the vans are to be stuffed. An of- 


fering is made based on the following formula: 


Booking _ ( Cargo- ) 


= ( Cargo ) Forecasted ) 
Oitering on-hand 


Booked aa Cargo Receipts 


The CSS computes the previous three weeks average weekly 
receipts and uses this figure as a basis for forecasting the 
Volume of cargo arriving at the facility between the date of 
Meeering and the stuffing of containers. ‘Figure 3 illustrates 
the present method. 

Figure 4:depicts the situation where a particular POD 
has 100 measurement tons (MT's) of cargo on the floor on 
Julian calendar day 70. To keep the computations simple, 
BeveMi*s of cargo will represent the average weekly receipts 
for this POD for the past three weeks (15 days, since only 
weekdays are included) giving average daily receipts of 50 
Mes. Container space has already been booked aboard three 
Meeocls during the 21- day period under consideration. 

The cutoff dates have been made for days 77, 79 and 84 for 
oe SO) and 120 MI's respectively. Using the past three 


weeks average receipts as the basis of the forecast, there 
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will be 360 MT's on-hand on day 77. The two-day period 
where no cargo accumulation occurs represents a Saturday 
and a Sunday. The five working days will accumulate 250 
MT's. Since space for 360 MIT's has been previously booked, 
all cargo on the floor for this POD will be stuffed. Two 
m@eomlacer there will be 100 MI's on the floor but 20 MT's 
mmerenot be stuffed. On day 84 only 120 of the 170 MT's 
will be stuffed. By day 91, the day for which the forecast 
meade, the forecast says there will be 300 MT's on the 
floor for this particular POD. The number of vans requested 
will depend on the average eube VeilnZae hone hr ene POD. 

If the POD's average is 50 MI's then an offering for six 
vans will be made on day 70 to be available for stuffing 


eee day 91, 


See CANCELLATION PROCEDURE 

Cancellation of space booked aboard a carrier must be 
made no later than "a reasonable length of time" prior to 
Mae cutotft date according to MSC's Container Agreement and 
meee ouide, RGS. If cancellation cannot be made within this 
meme and if the carrier cannot utilize this space, then the 
government may be charged for the space even though it was 
moteused. The present policy at WAMIMC is to cancel no 


later miieerve days Dr1IOr £o (CULOTL , 


Poe OPERATIONAL RESTRICTIONS AND REQUIREMENTS 
Mciomlamne Noted that there are other operational re- 
Sertetions and requirements that affect stuffing and booking 


procedures, although a thorough understanding of them is not 
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Pemeemtial to this work. Briefly they are: 1) a first-in- 
fet -Out (FIFO) cargo stuffing procedure requiring a mini- 
Mum average FIFO performance of 80% for each POD, 2) shipment 
Meemaraties that are exceptions to the FIFO procedure, 3) a 
Tequirement for the timely handling of certain goods such 

as household effects, 4) a requirement for a monthly cube 
MetiizZation average of 75% with no less than 50% cube utili- 
matron to ensure that container utilization 1S more economi- 
feeethan break bulk shipping, 5) a restriction on cargo 
mixing such as Military Assistance Program (MAP) cargo not 
being permitted to be mixed with any other type of cargo, 
Seume requirement of contracting the low cost carrier if 
Emere 15 more than one that can provide the required ser- 
meee and 7) a ceiling of 4,000 measurement tons on the floor 


at any one time. 
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Pi oiieobuLCIIVES AND PROCEDURES 


meee olLUDY OBJECTIVES 

The essence of this study is to observe the flow of cargo 
into the Container Freight Division at MOTBA and attempt to 
forecast by statistical and analytical techniques the expected 
Saneo inputs at some future point in time. The forecast is 
of the volume of containerizable cargo passing through the 
MrUeadestined for a particular POD. Break bulk cargo volume 
and source stuffed vans are not considered in this analysis. 

The necessity for such a forecast evolves from the fact 
that the Military Sealift Command must book container space 
aboard vessels as much as three weeks in advance of lift. 
Because of the demand for container shipping in recent years, 
commercial carriers are able to make this demand on their 
Customers. It remains to be seen whether a shift in the sup- 
my and demand factors will alter this requirement to any 
Sienificant extent. 

Although forecasting and forecasting models are generally 
associated with inventory and production control, the same 
Peiieral Ppranciples apply to the situation at MOTBA. Es- 
sentially, the situation is one of demand -- the demand by 
overseas commands for material support. In this respect the 
forecasting problem can then be broken down into components 
of that demand. Examples of these components are average 
demand, trend effects, seasonal effects, and noise or random 
EireetsS. by recognizing these components one can then con- 


SembuGErstatistical models to déal with them. 
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Accurate forecasting may have an influence on some of 
the performance variables that the container booking system 
generates. These are frequency of cancellations due to 
overbooking, age of cargo at stuff, cube utilization, and 
total volume of cargo on the floor. 

Mimeorecastine, the past is projected into the future. 

A good forecasting model should provide information feedback 
accurately and quickly so that the manager can make the 
Necessary and appropriate adjustments. It should provide a 
means of reflecting sudden shifts caused by policy changes 
within the system or incontrottaple exogenous factors which 
may also influence the system. A forecast that is based on 
out-dated information will not be beneficial to a manager 
who must book container space three weeks in advance. 

A forecast is only as good as the data it uses. [In this 
thesis it is assumed that the historical cargo flow data 
‘used to evaluate the forecasting methods is accurate. Data 
totals obtained from the computer tapes were verified against 


the totals that were manually computed at MOTBA. 


mee 6=PROCEDURES 

As stated previously the system variable that this thesis 
@earesses 1S the volume of containerizable cargo for a POD 
flowing into the Container Freight Division at MOTBA. Al- 
m@ough there are weight restrictions to consider, these re- 
merictions have not been significant factors in the booking 


mreocedure. 
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The data used throughout this thesis was taken from 
meem@al operations at the Container Stuffing Station. The 
data covers an 18-month period dating from 1 January 1973 
to 30 June 1974. This period generated an enormous volume 
of data which reflected the magnitude of the container 
stuffing operations at MOTBA. 

In order to handle this volume of information the ana- 
lysis was done in four separate phases. These were: 

1) a distribution audit in order to identify the vari- 
ables needed in the analysis, 

Eee che grouping of data in a more aeabile format, 

3) atest for serial correlation in the data, and 

4) the application of various forecasting techniques 
and their evaluation. 

MimperiLorming the distribution audit some of the com- 
puter programs from the McCarthy and Carter thesis were 
Meriazed. Only slight modifications to these programs were 
needed to manipulate the enormous volume of data generated 
by the eighteen months of operations. The audit identified 
the number of PW ioiemmcmanca tieevolume of these shipments. 
From this information one could identify major Ports of 
Weparkation, thus narrowing the scope of the problem. 

immorcaer tO arrange the data into a more usable format 
a number of short computer programs were written. These 
programs grouped the data from daily EOmweemry totals in 
the aggregate as well as POD levels. 

Since this work dealt with time-series analysis it was 


measonable to test the data for Serial correlation in order 
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to limit the range of potential forecasting methods. A 
Semeucer program was utilized to calculate the statistical 
information needed to evaluate the various PODs. 

Several computer programs were then used to apply various 
forecasting techniques to the data. The simplest techniques 
were applied first. In some cases these initial techniques 
were expanded to take into consideration the components of 
demand previously mentioned. The various forecasts were 
mem evaluated based on the difference between the fore- 
casted and actual volume. 

All the computer programs used in the data analysis were 
written in the FORTRAN IV computer language. These programs 
Were run on the IBM 360/67 computer at the W. R. Church 
momputer Center, Naval Postgraduate School, Monterey, 


@arifornia. 
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IV. DATA ANALYSIS 


A. DISTRIBUTION AUDIT 

Ban data was taken from a set of 19 computer tapes sup- 
plied by WAMTMC and was combined with existing data from 
previous research work to form a set of data covering 18 
months. The period covered was Quarters III and IV of FY 
"73 and all four quarters of FY '74. Computer programs 
Pemmenmused to EXtract the pertinent information from the 
tapes, combine it with existing Gace iden cee sehe |aco ie — 
gate information in alphabetical order by POD and consignee 
mia data file. 

Pomputer Program Two of the McCarthy and Carter thesis 
was used to delineate the data for further analysis. This 
peogran determined the number of shipments, the total volume, 
the average volume, the standard deviation of the volume, 
the total weight, the standard deviation of the weight, the 
mmeragpe density and the standard deviation of the density 
for each POD. The program also grouped each POD's shipments 
mreo intervals of volume, weight and density. 

iitemmicext procedure ECntailed the G€xXamination of the gross 
Peewly totals of shipments and volume of Cargo arriving at 
the container stuffing station during the 18-month period. 


The following information was obtained: 


SHIPMENTS 

Mearters [ITI & IV, FY '73 111,463 shipments 
rer ' 7 4 124,592 shipments 
18-month Totals Coo, Uso sii pments 
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VOLUME 


Miwarcers [II G IV, FY '73 Seovoago! Cubbie Leet 
my '74 O07 0),724 cubie feet 
18-month Totals 18,469,655 cubic feet 


There was obviously a tremendous amount of activity at 
MOTBA during the second half of FY '73. Over 47% of the 
meeal shipments and 45% of the total volume for the 18-month 
period occurred during the first six months of the data set. 

It was concluded that residual Vietnam conflict effects 
were responsible for generating this high volume of activity 
during the first six months of the 18-month period and that 
meamcid not represent the present and the expected future 
flow of cargo at MOTBA. Therefore, it was decided to use 
Omry the data for FY '74. 

The audit was run a second time using only FY '74 data. 
fieine that period 157 PODs and 3,356 consignees were 
identified as having received shipments during this period. 
These PODs and consignees accounted for 124,592 shipments 
fern a total volume of 10,070,724 cubic feet as previously 
Stated. This was an average weekly total of 2,396 shipments 
with an average weekly volume of 193,668 cubic feet (or 
4,842 measurement tons). Total weight for the period was 
[o5213,3335 pounds. 

ifereviewing the yearly gross totals for the 157 PODs, 
meewas obvious that 10% of these dominated the activity at 
MOTBA. Sixteen PODs were identified as being the most ac- 
Eive In terms of volume of cargo arriving at the container 


Stuffing station during the period. The cutoff point was 
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meettrarily chosen at 100,000 cubic feet of volume per year. 
As a result the least active of the 16 major PODs averaged 

a volume input of 2,041.35 cubic feet per week (or approxi- 
mately S51 measurement tons). The number of shipments for 
these major ports totalled 115,732 or 92.89% of all ship- 
Memos £or the period. Their volume total of 9,370,693 cubic 
feet was 93.05% of the total volume. 

Table I, Appendix A, shows the weekly gross totals of 
Shipments, volume and weight arriving at the CSS in the second 
meer Of FY '73. The totals for FY '74 appear in Table II. 
Table III, Appendix A, shows the data summary of the dis- 
fee DuUt1ion audit. Tabular data for the 16 major PODs is 
Mieesented in Table IV. Table V lists the activities of the 
major PODs by shipments, volume and weight. 

ine wcomputer output of the intervals of weight into 
which each of the shipments fell provided a look at the por- 
tion of the release unit shipments arriving into the system. 
It the number of release units were significantly large, it 
may be possible to forecast with greater accuracy the volume 
Saecargo expected to arrive during a specified period. 

There were 1,900 shipments of ten thousand pounds or 
meme during FY '74 using only data generated by the 16 major 
PODs. These shipments should have been classified as re- 
lease unit shipments and it is assumed that they were. They 
represented 1.64% of the shipment totals for the major ports. 

It was concluded that the amount of release unit ship- 
Ments arriving at the container stuffing station was insignifi- 
cant in view of these figures. Furthermore, operations at 
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MOTBA indicate that the expected arrival date of RU cargo 

was an unreliable item. Therefore, RU cargo was considered 
part of the uncontrolled cargo flow and grouped with LRU 

data for the purpose of this work. Table VI, Appendix A, 
Shows the number of RU shipments classified because of weight 
for each of the major PODs. 

It should be noted that no attempt was made to determine 
Sieenumber of shipments classified as RU for reasons other 
than weight (1.e., classified material). It was assumed 
that the volume of this cargo was of insignificant proportion 
when compared to the total reensly volume. 

In summary, the distribution audit provided a means of 
taking raw data and arranging it into a meaningful and use- 
miiecormat. It permitted the identification of the most 
active ports of debarkation which then became the focal 
meint £Oor further analysis. The grouping of the data into 
weekly intervals by POD indicated the tremendous volume 
imeictuations that occurred throughout the period and the ap- 
Parent complexities that arise when attempting to forecast 
future volume flow. 

ipmeomeni1s point, only the simple computer manipulation 
of data had been accomplished through the diserfimeton aUiCueee 
With this phase completed, there were still many questions 
to be answered. What type of data had been generated? Was 
there any correlation in the volume flow from week to week? 

The above questions are answered in the following section 
foren Degins the statistical analysis of the data. The ans- 


wers then provide statistical justification for applying 
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Semedit Lorecasting techniques which are covered in section 


Weer this chapter. 


Pee oe RIAL CORRELATION TEST 

ime flow of uncontrolled cargo volume into the Container 
Poemont Division at MOTBA is time-series data. It is not 
iieannon Lor future observations in real time-series data to 
memccpendent upon current observations. This dependence may 
be measured by analysis of the consecutive disturbances around 
a Simple least-square-regression line. If these disturbances 
are found to be correlated from one observation to the next, 
Maaecnus not random, serial correlation is said to be pre- 
mae it the presence of serial correlation can be properly 
miimer.1ed it Can result in a more efficient estimation 
Peoeecss In comparison to simple least-squares estimators. 
euened process would give more weight to recent data for 
meeacime predictions of the values of future observations. 

Johnston [Ref. 9] lists three main consequences that 
result from applying straightforward least-squares formulas 
mmpect ly to observations that contain serially correlated 
Mesturbances. They are: 

Meeeline Sampling variances Of the intercept and the slope 
Seertticient may be unnecessarily large in comparison to those 
Pecained by other methods. 

Dye the usual least-squares formulas for the sampling 
Pirainces sor the regression coefficients are no longer valid 
Pieces a serlous underestimate of the variance is likely to 


SeCcur, 
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Sieetic sanplana variances of the prediction will be 
Mnnecessarily large. 

1. The Durbin-Watson d-Statistic 

When making forecasts using time-series data it 

would be a serious error to assume serial independence of 
the disturbance term. In order to avoid such an error, the 
Durbin-Watson d-Statistic [Ref. 6] provided a suitable test 
momeene presence of serially correlated disturbances. The 


Durbin-Watson d-Statistic is defined by 


nh 2 
Lo (Z, - Z._4) 


qe 
5 Line 
t=1 (2+ 


where Z,(t=1,2,...,n) denotes the residuals from a fitted- 
least-squares regression. To determine positive serial 
correlation the observed value of d is compared against 


lower and upper bounds of d, and d,, that are tabulated for 


L U 
various values of n (the number of observations) and k (the 
mamber Of EXplanatory variables}. The d-statistic for nega- 
mame serial correlation is calculated by taking 4 - d. 

This figure is also compared with dy and diy: SIIMeES ele Saves 
Of the serial correlation was unknown, two-sided tests were 
momdueted by combining single-tail tests. The test has four 
mossible outcomes: (1) positive serial correlation of the 
ae 2) Negative serial correlation of the data; (3) in- 


@oneclusive evidence of serial correlation with further 


Observations being ideally required; or (4) no presence of 
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memeai correlation (which implies that the disturbances are 
random). 

The calculations were made using three different 
intervals. The first interval used weekly data resulting 
memozeobservations. The second interval used biweekly 
Pemals beginning with the observation for the first two 
weeks of FY '74. The total number of observations in this 
set was 26. Biweekly totals beginning with the observation 
eeete Second and third weeks of FY '74 (1.e., ignoring the 
first week's volume total) was used for the third interval. 
Thus there were only 25 observations resulting from this 
meeuping scheme. 

The following example may help to clarify the above 
procedure. Given MidttiriomVvCcial) sVeMuMcmtipiuc for POD, BAS, 
momecnhe first seven weeks of FY '74 is 10, 20, 30, 40, 50, 
mumama. 7) cubic feet respectively, the first seven observa- 
meems tor the weekly interval data set correspond according- 
mameeine first three observations for the first biweekly 
data set will be 30, 70 and 110 cubic feet. However, the 
ijerst three observations for the second biweekly data set 
wall be 50, 90 and 130 cubic feet. 

The d-statistic test was implemented through Com- 
mieereerogram One. In addition to displaying the results 
feene calculations, the program displays the meaning of 
DitMmecmitlistte (1,¢., positive, negative, inconclusive or 
iii computer Output One at the end of this thesis is a 


Sample ser the printout for Computer Program One. 
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meeeSUlts OF ocrial Correlation Test 

While the Durbin-Watson d-Statistic was calculated 
for all 157 PODs, the d-statistic for: the 16 most active 
PODs was the focal point for the analysis. 

The results of the test for the weekly interval data 
set having 52 observations for each POD indicated positive 
mental correlation in the weekly volume inputs for nine PODs. 
The test was inclusive for four PODs and indicated no serial 
@emrelation for three others. There were no PODs with nega- 
Meversecrial correlation. In fact, there was no indication 
of negative serial correlation under any of the circumstances 
tested for all 157 PODs. The following table shows the test 
results using only the 16 major PODs. The test was conducted 
foetekevel of significance of 0.05. 


Week per Week nee No 


Interval Started Positive Negative Incon Correlation 


1 1 9 0 4 5 
Z 1 5 0 5 6 
Z 2 7 0 5 4 


manlesylingives the d-Statistzre results for each POD 
fmemechne three intervals tested. 

Crvenetne results of the d-statistic tests it was 
mmemeaecided to test for the aggregate presence of serial 
memerelation in the data for the 16 major PODs. The null 


hypothesis for the test was 
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where Pyr-s+ Pye Someday rtcopmesents ene serial ecorrela- 
meemecoctticients for the 16 major PODs. The alternate 
hypothesis was 


Hy: not all o's = Q. 


iars hypoemesis 15, in essence, testing the plausi- 
fmety that the individual cases indicating positive serial 
Semel atlon were Type I errors, i.e., a rejection of the 


null hypothesis, p. = 0, when it was true. This was formu- 


i 
lated as a Bernoulli process with the probability of success 
(the indication of positive serial correlation) equal to the 
momel of significance for the d-statistic test (0.05). 

Since a Bernoulli process was involved, the binomial 
Oestribution function was applied. Letting the probability 
Memmetieccess, p, equal 0.05; the number of trials, n, equal 


to 16; and the specific number of successes, r, equal to 5, 


feema 9 respectively, the decision rule was 
Accept H, if P(R>r|n,p) > 0.05 
Reject H, if P(R>r|n,p) < 0.05 


where R equals the unknown number of successes. 
verme Cumulative probability distribution tables the 


mollowing results were obtained: 


P(R>5|16,0.05) = 0.0009 


0.0001 


P(Re7 | Voeh05 ) 


P(R>9]16,0.05) < 0.0001. 


Based on these results the null hypothesis was re- 


meetea and it was assumed that significant serial correlation 
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Surely €Xisted somewhere in the data. Therefore, there was 
justification for using methods other than ordinary-least- 
Sauares regression to forecast the volume of cargo arriving 


mune Container Freight Division at MOTBA. 


C. FORECASTING METHODS 

Moree results of the Durbin-Watson d-Statistic indicated 
mmempresence of serial correlation and thus the serial de- 
mema@ence Of the observations. It showed that the distur- 
Dances around an ordinary-least-squares regression line 
Pomme not fall equally on each side of the line. As the 
minturDances tended to fall more on one side of the line 
mmemeche other, this implied that the most recent distur- 
Mimeeswes the best predictor of the subsequent distrubance. 
mimeo y amplication, the best predictor of the next obser- 
mreemom 15 Che most recent observation or observations. This 
Seeplied justification for the application of forecasting 
meenwods Other than ordinary-least-squares to place more 
EMpiasis on recent data rather than the minimization of the 
eared residual term. While there were a number of methods 
from which to choose, the methods used in this work were 
mmesem because of their relevency to the problem and for 
miedr ease of calculation. 

This section looks at: the moving average method that 
is presently being used at MOTBA; the exponentially weighted 
average method; the exponentially weighted average method 
Meieadyusement for trend; and an adaptive response rate 
mcnmrque applied to the exponentially weighted average 


method. 
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These methods were evaluated using mean absolute devia- 
tion (MAD) and mean deviation squared (MDS) as the criteria 
for goodness. By computing the mean of the deviations it 
was possible to compare the various forecasting methods when 
they resulted in differing numbers of forecasts generated 
during the period. 

In addition, MAD and MDS were computed for one-, two- 
miceethree-week time horizons. For the one-week time horizon 
the forecast was subtracted from the actual input volume on 
a week by week basis. For the two-(three-)week time horizon 
the latest forecast was nner by two (three) and sub- 
fjeacted from the actual input volume for the two (three) 
weeks covered by the forecast. This was done on a week by 
week basis. 

ives tOollowing example may help to clarify time horizon 
mebenlations. Say that the forecast for week 20 is 1,000 
Meercmrcet. If a three-week time horizon 1s being used, the 
forecast will be 3,000 cubic feet for weeks 20, 21 and 22. 
memeie adctual input volumes turn out to be 900; 1,100 and 
meee cubic feet respectively, the absolute deviation for 
Maer torecast is then 200 cubic feet (or 900 + 1,000 + 1,200 
- 3,000). Week 21 will result in a new forecast for weeks 
mez e and 25. Ihe forecast for week 21 is multiplied by 
three and subtracted from the actual input volumes for weeks 
meeez2 and 25 to determine the absolute deviation for that 
Paeemeuvdr forecast. thus, for each forecast there 1s an 


absolute deviation that is computed. 
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What, then, was the purpose for using time horizons? 
Siremeteason was the requirement for booking container space 
three weeks in advance and thus projecting a one-week fore- 
cast over a period greater than one week. One-, two- and 
three-week horizons were used to make comparisons of any 
effects that may be caused by changes in the booking re- 
quirements. It was also anticipated that an increase in 
Piemctame horizon may have some effect on the results of the 
various forecasting methods due to a dampening effect on 
G@etaewith large fluctuations. This dampening effect is ex- 
plained later. | 

It should also be noted that comparisons among the var- 
ious forecasting methods were only made within a given time 
Meordzon. 

Tables VIII and IX, Appendix A, show the MADs and MDSs 
generated by all of the forecasting methods used in this 
analysis. 

1. Moving Average Method 

As previously explained, the CSS at MOTBA computes 
the previous three weeks average receipts in order to fore- 
cast the input volume that will be on hand in 21 days. The 
mercrinegs for bookings that are made depend on three factors: 
Mmetne volume of cargo on hand; 2) the volume of cargo 
booked; and 3) the forecasted volume. 

The three-week moving average is nothing more than 
feline ar average. This average, multiplied by three 1s then 


Heed to £orecast the volume receipts for the next three 
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weeks. In this method equal weight is given to each of the 
past three weekly volume totals. 
In order to evaluate the moving average method the 
MADs and MDSs were calculated for one-week to ten-week 
moving averages for the 16 PODs. In addition, the aggre- 
gate MAD and MDS for each POD were calculated for one- to 
ten-week moving averages. In other words, a specific moving 
average (from one to ten-weeks) was applied to the volume 
mapucs for each POD. The MAD and MDS for each POD were 
@atculated and then summed to get the aggregate. The same 
@atculations were made for two- and three-week time horizons. 
Computer Program Two was utilized to obtain the 
calculations for the analysis of the moving average method. 
Computer Output Two shows a portion of the output that was 
generated by the program. 
a. Data Analysis of the Moving Average Method 
In 13 of 16 cases a moving average other than 
Mmfemicnree-week Moving average provided a better forecast 
when MAD was the criterion of goodness. When using MDS as 
Mmewecriterion, the result was 15 of 16 cases. The time in- 
fPamalse that provided the best forecast will be called the 
"optimal moving average." 
These evaluations were made using a one-week 
Mame MOrlizOn. Similar results were obtained using two- and 
Mittecowechectime horizons except that the optimal time inter- 
val for the moving average was not the same in each case. 
It was noted that in some cases the selection of 


Piomougganess —eriterion determined the optimal interval for 
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Biesmoving average. This was the case for six PODs using 
a one-week time horizon and was true to a greater or lesser 
extent for the two- and three-week time horizons. 

The decrease in the MAD when a moving average 
other than a three-week moving average was optimal, varied 
meomea,o0 tO 18.5%. Similar reductions in MAD and MDS oc- 
curred when the other time horizons were used. 

The optimal time interval for a moving average 
that applied to all PODs was seven, ten and ten weeks for 
Smee twO- and three-week time horizons, respectively, with 
MAD as the criterion. The same time intervals held true 
meenerDS as the criterion. 

It was noted that the number of forecasts de- 
Greased as the time interval for the moving averages in- 
meericed, For example, in a three-week moving average, the 
first three weeks are averaged to make a forecast for week 
four. Thus, there were only 49 forecasts made for the 
meamroa. For a ten-week moving average there were only 42 
memecastS. The use of two- and three-week time horizons 
mesulted in even fewer forecasts. Comparison among the 
various moving averages for a specific time horizon was 
possible since mean absolute deviation and mean deviation 
squared were used as the criteria for goodness. 

Dee sine winter valesor the mMovin® average i1n- 
creased, the greater was the chance for large fluctuations 
mimcewe; data to be dampened by the initial forecast. For 
example, the volume inputs for the first six weeks of FY 


peomeoemewewere 2,611; 4,307; 1,079; 3,360; 6,765 and 
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@yreo Cubic feet, respectively. These fluctuations were 
dampened when moving averages greater than three weeks were 
applied to the data (with a five-week average being opti- 
mal). This dampening effect stressed the fact that a three- 
week moving average placed more emphasis on recent data. 
memeer time intervals placed less emphasis on recent data 
and considered older data as well. 
PeeWee Soenotcdstrat the MAD and MDS decreased 
as the time horizon iMmegeascd as (ne reasonewas. that the 
two- and three-week time horizons also tended to cancel 
the volume fluctuations. A look, again, at Figure 1 shows 
fmemvery large volume fluctuations that occur in the POD 
data and how it is possible for the dampening effect to take 
place. 
fee Exponentially Weighted Average Method 

iimeeolweor thesresults Of the Serial correlation 
mest, a Viable alternative to the moving average method is 
the exponentially weighted average (EWA). In its simplest 
form an EWA is based on a period by period adjustment of the 
most recent forecasted average by adding or subtracting a 
Braction of the difference between the actual volume in the 
@emrent period and the latest forecasted average volume. 
memditters from the moving average in that the weight given 
the most recent observation can be varied while also giving 


Sonsideration to all past data. EWA is represented by 


Eee ee ey) ee 
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where Be Hoephewtorecast average, V¥. 15 the actual input 


iE 
volume and a is the exponential smoothing constant which 


must be between zero and one. The forecast for the next 
* 


period F.., 


can be taken directly from the computed value 
of F,. This can be justified because trends and seasonal 
adjustments are not accounted for in the model [Ref. 7]. 

The exponentially weighted average method is also 
known as the exponential smoothing method. Both names are 
used in this thesis. 

The calculations for the EWA method were implemented 
by Computer Program Three. The MAD and MDS for each POD 
mene calculated for a-levels between -0.10 and +1.00 in in- 
crements of 0.01. Yn addition the MAD and MDS were computed 
Mmoeeo-levels that were applied to all 16 PODs. In othtér 
words a specific a-level was applied to data for each POD. 
The MAD and MDS for each POD were calculated and then summed 
Pomobtain an aggregate total. The a-levels which generated 
the smallest MAD and MDS are referred to as the "optimal 
momevels for aggregate data." It should be noted that the 
optimal a-levels for aggregate data may differ depending on 
fmm criterion of goodness. 

Pmeciied COMpULcCrEESubDrourine, ocUBROUIINE OSPLOT, 
available at the computer center, was used to plot the 
a-levels versus the MAD for each POD. The plot displayed 
how the MAD reacted to various levels of a. Thus it was 
messtolec te Visually determine the optimal a-level for 


aoe race data with MAD as the criterion of goodness. 
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Computer Output Three-Two shows the MADs and MDSs 
Calculated for various a-levels. Computer Output Three-One 
1s a sample of the OSPLOT for four of the PODs. 

a. Data Analysis of EWA Method 

With MAD as the criterion for goodness, the EWA 
method resulted in a decrease in MAD in 14 of 16 cases when 
compared to the three-week moving average. The decrease 
ranged from 4.9% for RJ3 to 12.4% for XJ1. However, when 
EWA was compared to the optimal moving averages for the 
Mea dual PODs it produced a smaller MAD in only 4 of 16 
cases. The a-levels varied from 0.04 to 0.60 for the 16 
PODs. 

iiewopeiicme -lMevelerorsaporegate MAD data was 
0.36. The aggregate MAD for the EWA method was 59,387 
Smee teet as compared to 60,172 cubic feet for the three- 
week moving average and 59,209 cubic feet for the optimal 
Meving average. 

The preceding analysis was made from calcula- 
tions using a one-week time horizon. Similar results were 
Sptained using two- and three-week horizons. It is interest- 
meoeto note that there were changes in the a-levels for most 
memeoeaepending on the time horizon. One POD (UM4) produced 
megative a-levels (-0.02 and -0.01) for two- and three- 
week time horizons, respectively. The optimal a-levels for 
aggregate data were 0.35 for the two-week horizon and 0.33 
fOr the three-week horizon. AS was the case with the moving 
average method, the aggregate MADs decreased as the time 


Pearezolvsmincreased due to the dampening of the fluctuations. 
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With MDS as the criterion for goodness the EWA 
method resulted in a decrease in MDS a 12 of 16 cases when 
compared to the three-week moving averages. When compared 
to the optimal moving average, there was a decrease in only 
5 of 16 cases. In each case the optimal a-level for MDS 
was different from the optimal level for MAD; although the 
optimal a-level for aggregate data was approximately equal 
for the MDS (0.37) and the MAD (0.36). The aggregate MDS 
for the EWA method was less than the MDS for the three-week 
Ferenc amerage DU greater than the MDS for the optimal 
moving average. 

A one-week time horizon was used in the above 
analysis for MDS. As was the case with the MAD analysis 
of individual PODs, two- and three-week time horizons 
penerated similar results. However, for a two-week time 
horizon of aggregate data, EWA generated a slightly smaller 
MDS than that generated by the optimal moving average. 

This was not the case for the three-week time horizon, 
however. 

The way in which the EWA method was applied re- 
feirted in no forecast deviations being computed until the 
mimercdeweek. this resulted in some smoothing of the data but 
not as much smoothing as occurred with moving averages with 
Mmicervals greater than three weeks. When applying the EWA 
Mmerrvod to data with large fluctuations at the beginning of 
miemyvedar, there was a lag in the forecast until these fluc- 
tuations were dampened over time. This was the apparent 


reason for the seven-week moving average being a better 
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menecaster than the EWA for a one-week time horizon and a 
ten-week moving average being better for a two- and three- 
Week Cime horizon. 

In summary, the EWA method provided a means of 
emphasizing past data as well as placing some significance 
on the most recent data by the use of the smoothing con- 
=eae, the EWA method provided a better forecast than the 
manree-week moving average; however, the decrease in the MAD 
Pexeaperegate data was only 2.2%. For individual PODs this 
decrease was as high as 18%, but for other PODs there was 


an increase in the MAD when EWA was used. 


3. Exponentially Weighted Average Method with Trend 
Adjustment 


\Vipicm themmovineeaverage or the EWA took into 
consideration any adjustment for trends that may exist in 
eiemaata. Also, the slopes of the least-squares-regression 
line for some PODs used in the d-statistic calculations gave 
fmerndication of the possibility of trend effects in the 
data. Therefore, the EWA was modified to take into account 
Sg@emadajustment for trend. 

The difference between the volume forecast averages 
from week to week is the apparent trend. Exponential smooth- 
mececan be applied to this trend just as it was applied to 


the forecast average. This trend adjustment is represented 


by 
ee eer Clo) T, 4. 


Byscombinine this trend adjustment with the new forecast 


avetace, the forecast for the next period is then 
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= = 1 = 
Piper a) Eee esa 
or 
= hee TF 
Poe tO OE 
where F* 1s themtorecast for the exponentially weighted 


et | 
average with trend adjustment (EWA-TA) [Ref. 7]. 


Computer Program Four was used to make the statisti- 
@meecalculations for the EWA-TA. Except for the trend ad- 
Ererment , the computational procedures were the same as those 
for the EWA method. Once again, SUBROUTINE OSPLOT was used 
memolot a-levels versus MAD for each POD. 

a. Data Analysis of EWA-TA 

Using a one-week time horizon and MAD as the 
mimtetTi0On, a decrease in MAD resulted in only 3 of 16 cases 
when comparing EWA-TA to EWA. Only in the case of RG1l was 
fiaere any significant change (an 8% decrease). However, 
mse’ OD was an exception to the rule of fluctuating volume 
inputs. Its weekly volume inputs toward the last half of 
maceyear had remained within a relatively small range. In 
Sy one case did EWA-TA provide a better forecast than any 
fete methods previously discussed. However, the improve- 
Meme was minimal. In half the cases the difference in MADs 
between EWA-TA and EWA was less than 50 cubic feet. Also, 
mg@ema-=level dropped significantly in all 16 cases. The op- 
timal a-level for aggregate data was 0.16; however, the MAD 


Mes greater by 1,904 cubic feet. 
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With MDS as the criterion for goodness the re- 
Sees were similar. Improvement on the EWA method occurred 
in only two cases. In most cases the a-levels were differ- 
ent from those with MAD; however, the optimal a-level for 
aggregate data changed by only 0.01 to 0.17. The aggregate 
[@isewas the largest of all the methods applied thus far. 
This was true for all time horizons. 

Very few changes occurred when two- and three- 
week time horizons were applied. Once again, UM4 produced 
negative a-levels. The optimum a-levels for aggregate data 
did not change significantly. The aggregate MAD and MDS 
decreased in the same proportion as the aggregate totals for 
maceother methods. 

In summary, the exponentially weighted average 
method with trend adjustments did not provide an improved 
forecast over the moving average method or the exponentially 
weighted method without trend adjustments. Apparently there 
memo derinite trend in the volume inputs for the PODs. It 
memaoubtful that the method can be effective when applied 


memcata covering only a one year period. 


4. Exponentially Weighted Average with an Adaptive 
Response Rate 


Because of the large volume fluctuations in the data, 
the application of an adaptive response rate (ARR) provided 
a plausible alternative to a simple exponentially weighted 
average. The ARR provides an automatic means of detecting 
Pacaen volume changes through the use of a tracking signal. 


Mmimceeana beach [Ref. 18] define the tracking signal as 
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Smoothed Error 


Tracking Signal = Smoothed Absolute Error 


where the error is the difference between the forecasted 
volume and the actual volume inputs in each period. The 
error was smoothed by choosing a value for gamma, y, between 


=imeand +1 such that 


New Smoothed Error Give old snooemed error + yy error. 


New Smoothed Absolute a 


Peau (1l-y) old smoothed absolute error 


+ y absolute error. 


The tracking signal was then applied to the EWA methods with- 
out trend adjustment by letting the smoothing constant equal 
mftemabsolute value of the tracking signal. Thus, it ap- 
peared that ARR would provide a rapid response to large 
volume changes whereas the constant model would lag behind 
these sudden changes. 

The adaptive response rate was also applied to the 
exponentially weighted average with trend adjustment. It 
resulted in very large deviations because of the conflict 
meeween the trend adjustment and the response rate. They 
were apparently working against each other; and therefore, 
mies method will be ignored in this paper. 

MADS and MDSs for the 16 major PODs were calculated 
for y-levels between -1.00 and +1.00 in increments of 0.01. 
Optimal y-levels for aggregate data were also computed. 


These calculations were implemented by Computer Program 
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Fave. Plots of y-levels versus MAD were made through the 
use of the special subroutine 
a. Data Analysis of EWA with ARR 

With a one-week time horizon and MAD as the 
Gre rion, the ARR resulted in the smallest MAD in 6 of 16 
cases when compared to all other methods previously dis- 
Emesced. Ihe reduction in MAD was less than 5% in five of 
these cases. The optimal y-level for aggregate data was 
-0.05. The aggregate MAD at this level was smaller than 
the MADs for the three-week moving average and the EWA-TA; 
Wipeett Was larger than the MADs for the optimal moving 
average and the EWA method. 

A two-week time horizon produced similar re- 
sults on the individual POD as well as aggregate levei. 
im@evidiual POD results were similar with a three-week time 
horizon. For aggregate data only the optimal moving aver- 
age resulted in a smaller MAD. However, the reduction was 
Mess than 17, and for all practical purposes it was equal. 

With a one-week time horizon and MDS as the 
criterion, the ARR produced a lower MDS in 5 of 16 PODs 
rien compared to all other methods. Two- and three-week 
mnenmiortzons had similar results on the POD level. For 
aggregate data, the ARR was unable to generate an MDS smaller 
meomeany Of the other methods for the three time horizons. 

In summary, the ARR provided a more complicated 
approach to exponential smoothing than any of the other 
methods. On the individual POD level it provided a better 


forecast in some cases; yet, the improvement was minimal. 


On 





i@emcime horizon and the criterion for goodness were signif- 
icant factors in determining whether the ARR generated the 
better or best forecast. 

While the purpose of the ARR was to respond more 
rapidly to large fluctuations in the data, it would be more 
mepnopriate to state that it was designed to respond to 
meee Fluctuations which resulted in a shift in the data. 
While there were large fluctuations in the data there were 
Eaomapparent shifts in the data. Plots of the weekly volume 
inputs indicated that there were large fluctuations in the 
fea, bUt that these fluctuations did not result in any 
Significant shifts in the data. 

5s. Summary of Forecasting Methods 

hour £orecastingm=techniques were applied to the 
meee 1) the moving average method; 2) the exponentially 
Memented average method; 3) the exponentially weighted 
menace method with trend adjustment; and 4) the exponen- 
tially weighted average method with an adaptive response 
mace. Past data was utilized to forecast future inputs in 
each of the methods. Each method provided a means of 
varying the amount Pec Npiac hom raectimMoOlmene most recent 
Mira, £Orecast average, or error. 

immerder to evaluate the forecasting methods, the 
mean absolute deviation and mean deviation squared were 
@abeulated and used as the criterion of goodness. The mean 
of the deviations was used since each method resulted in 


a different number of forecasts being made for each period. 
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The MADs and MDSs that were generated by the various 
forecasting methods indicated the uniqueness of the volume 
inputs for each POD. The method providing the best fore- 
cast for a POD depended upon the criterion of goodness and 
the time horizon. 

Tables VIII, IX, X and XI, Appendix A, are a pre- 
sentation of the data that has been previously discussed. 
Table VIII gives the MAD of the various methods for the 16 
major PODs in one-, two- and three-week time horizons. 

Table IX gives the MDS in the same format. Tables X and XI 
present the MAD and MDS for aggregate totals respectively. 

From the data in these tables it can be seen that 
lower MADs and MDSs may result from increasing the time 
horizon. It has been suggested that the lower deviations 
were a mathematical result of the two- and three-week time 
horizons averaging the volume fluctuations. The reduction 
Someimued for each increase in the time horizon in 8 of 16 
PODs. In only one case did a decrease not result from an 
mmerease in the time horizon from one to two or three weeks. 
mamother words, in 15 of 16 cases there was a reduction in 
MAD when using a two- or three-week horizon. 

The moving average method provided a means of damp- 
ening the initial large volume fluctuations in the data. 
isaree volume fluctuations during the first quarter of FY 
'74 could be dampened by using a 10-week moving average. 
mms, the first deviation was not computed until the eleventh 


week. 
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The other methods that were applied to the data be- 
gan computing deviations as early as the third week. With 
maeee VOlume fluctuations at the beginning of the period it 
was possible that there was a lag in the forecast which 
faeetally resulted in large deviations. The plots of the 
weekly inputs showed that 8 of 16 PODs had large volume 
mietuations at the beginning of the period. : 

An overall view of the forecasting methods indicated 
that no one particular method dominated the others in terms 
BE providing the best forecast on a POD or an aggregate 
level. The difference between the largest and smallest MAD 
Was less than 100 cubic feet for one POD (UL7). For another 
POD (UM4) three different methods were optimal depending on 
the time horizon. These and other factors demonstrated the 
mamcom Llow of cargo arriving at the CSS. 

Peeeeopecimal Forecasting Method Selection 

The results of the various forecasting techniques 
applied to the data provide a basis for the selection of an 
Serimal method. Through the application of these techniques 
the following alternatives were available for forecasting 
the weekly volume inputs of the Pops:?+ 

(1) a three-week moving average, 


(2) a moving average with an optimal time interval 


meplied to all PODs. 


Icee Tables X and XI, Appendix A. 
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(3) a moving average with the optimal interval for 
facie POD applied to its weekly inputs (i.e., the time in- 
tervals vary among the PODs), 

(4) an exponentially weighted average with the op- 
timal a-level applied to all PODs, 

(S) an exponentially weighted average with the op- 
femamea-level for each POD applied to its weekly inputs 
(i.€., the a-level varies among the PODs), 

(6) an exponentially weighted average with trend 
adjustment with the optimal a-level applied to all PODs, 

(7) an exponentially weighted average with trend 
adjustments with the optimal a-level for each POD applied 
to its weekly inputs, 

(8) an adaptive response rate with the optimal 
fmeeevel applied to all PODs, 

(9) an adaptive response rate with the optimal 
meee vel for each POD applied to its weekly inputs. 

The selection of the optimal method required an 
melysis of the methods used to calculate the goodness 
@edeerion and the application of Student's t-test in order 
to determine if the difference of the deviations among the 
memmeas were significantly different. 

The results of the one to ten-week moving average 
indicated that a time interval greater than three weeks 
generated smaller deviations than the three-week moving 
Mmeyeragc for 13 of 16 PODs. This ratio held true for both 


errcerions and for all time horizons. Thus, if a moving 
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average is to be applied to the data, a time interval greater 
than three weeks is preferred. 

There is one shortcoming in comparing the moving 
meemage method with the others. The first deviation for a 
moving average was not computed until the week following 
the time interval of the average while the deviations were 
computed beginning the third week for the other methods. 
Thus, the longer time intervals tended to have a greater 
imitial smoothing of the forecasts while there was a lag 
in the smoothing for the other methods. This fact should 
be kept in mind when comparing the aggregate totals in 
iaple X, Appendix A. Because of the manner in which the 
deviations were computed no comparisons were made between 
the optimal moving average and the other methods. 

iiieme xponemtctally welghted average otters an alter- 
Mmerve tO a moving average. Firstly, it did generate smaller 
MADs and MDSs than a three-week moving average. A simple 
femest using MAD as the criterion and a one-week time hori- 
memeindicated that alternatives (4) and (5) are significantly 
fimererent from alternative (1) at the 0.05 confidence level; 
Sereetaeretore, (4) and (5) are preferred to (1). At the 
Same confidence level, the t-test indicated that (4) and 
memeare different; and, therefore, alternative (5) is pre- 
moneda to (4). Not until a confidence level of 0.005 is 
Meedemean tt be said that alternatives (4) and (5) are not 
ecimiiitcantly different. 

The exponentially weighted average method with trend 


adjustment (alternatives (6) and (7)) generated aggregate 
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MADs and MDSs that were larger than those generated by the 
Preeewithout the trend adjustment. Based on this fact, al- 
Bemiatives (4) and (5) are preferred to alternatives (6) 
and (7). 

The exponentially weighted average with an adaptive 
response rate (alternatives (8) and (9)) provided another 
atermnative to the moving average method. For 6 of 16 PODs 
the ARR generated the smallest MAD for one- and two-week 
Bemenorizons. On the aggregate level alternatives (8) and 
(9) generated slightly smaller MADs than did alternatives 
maeand (5). A t-test indicated that alternatives (8) and 
(9) are different; and, therefore, (9) is preferred to (8). 
The test was conducted at the 0.05 confidence level. 

Wath alternatives (5) and’ (9) generating smaller 
MADS and MDSs than the three-week moving average, the ques- 
Mmeenewas which alternative is preferred? Again, a simple 
t-test at the 0.05 level of confidence was applied. At this 
level, the test indicated that alternatives (5) and (9) are 
Smeecaeniticantly different. With this fact in mind, al- 
ternative (5) would be preferred to (9) because of compu- 
tational ease.” 

Himlcomttewas possible Lo Use a sample t-test to ar- 
maverat a preferred method of forecasting, the results of 
the various methods have indicated the true randomness of 


m@mercargo flow into the Container Stuffing Station. No one 


‘Table XII, Appendix A, shows a sample of the t-test 
ealeulations. 
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Meneicular method dominates the others in its ability to 
meeurately forecast the weekly volume inputs for all of the 
major PODs. The satisfaction of arriving at a preferred 
method must come from the fact that this method does result 
in deviations that are smaller than those generated by a 
three-week moving average and is computationally easier and 
less sensitive to the data than another method that is com- 
parable in its forecasting accuracy. 

Additionally, the minimum MADs achieved using alter- 
PEisive (9) occurred at y-levels ranging from -1.00 to +0.49. 
Also, the data for most PODs possessed several maxima and 
minima of MAD as y-levels were varied. This revealed the 
Naeghly sensitive nature of the results of forecasts to the 
mmput data when alternative (9) was used. Should the char- 
@eteristics of these data change in future observations, 
Merers resulting from the continued use of previously opti- 
mal y-levels could become substantial. This was contrasted 
to the more narrow range of a-levels (from 0.02 to +0.77) 
momealternative (5) and the achievement of a single minimum 
MAD for most cases. 

icomoueet Gueput Ihree-One and Computer Output Five- 
Smienare representative of these relationships. 

fomerorecast Variance Estimation 

foments point, this thesis hasbeen concerned with 
mmeechoice of an optimal method of forecasting the volume 
inputs to the CSS, using mean absolute deviation and mean 
[mvrreyon squared as the criterion of goodness. The appli- 


Mieroneor any method chosen requires that the statistical 
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Variance of the forecast be known or estimated such that 
confidence in the forecast may be assessed. The knowledge 
of the variance can assist the booking of vans by providing 
a means of assessing the impacts and costs of conservative 
versus aggressive booking policies and the probable effects 
on average age of cargo, cancellations, cargo volume on the 
moor, cube utilization and other policies. Just as the 
expected volume was forecast by some method, so must the 
expected variance be forecasted or estimated. Two possible 
procedures for arriving at variances for the forecasts are 
suggested. 

The first possible method of estimating the variance 
and standard deviation of volume inputs for a POD could be 
bewproject forward historical information for one period 
to the next. If this method were chosen, the length of the 
period would dictate how much averaging out of fluctuations 
would be accomplished. Conversely, should the fluctuations 
be tied to some major policy decisions occurring during a 
Memwod, Variance estimates could be seriously in error if 
mmeneperiod chosen was too long. It is suggested that per- 
haps a period of six months to one year, providing 26 to 
S2 weekly deviation observations. Thus, the weekly volume 
mpi variance for a POD calculated for one period could be 
applied to all weekly forecasts for that POD in the ensuing 
period. 

Werlization of a smoothed square deviation computed 
Pinim@btamecously with each volume forecast is another possible 


method of estimating the variance. As the variance of each 
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forecast is related to the square of the deviation between 
Eiewobserved and forecasted volume for the period, this 
value could be smoothed from one forecasting period to the 
next using the simple exponentially weighted average tech- 
meaue. lhe estimated standard deviation for the next fore- 
Saaeecoulc be calculated as the square root of the smoothed 
Squared deviation. The value of the smoothing constant 
(alpha) for this calculation could be selected depending on 
the fement to be placed on the most recent data (usually 
merween 0.1 and 0.3 according to Brown [Ref. 3]). Intui- 
Mmmermy, the selection of a larger alpha could be tied to 
Meme aperessive booking policies by trying to detect recent 
trends while smaller values would provide more conservative 
/Meme-run averages. In any case, this technique provides 
continual updating of variance estimates. 

Neither of these methods was tested persuant to 
Mmmm thesis due to the limitation of having only one year's 
eed and the lack of any apparent measure of goodness that 
would indicate the relative optimality of one method over 
mies other. 

Whether or not either of these methods is used, 
some estimates of the variance associated with a forecast 
prediction is imperative to establish consistency between 
mresbooking policy and other policies concerned with the 
movement through and storage within the Container Stuffing 


mation. 
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Meee CRITIQUE OF PROCEDURES AND DATA ANALYSIS 

The analysis leading up to and including the forecasting 
Beeniniques have been presented in detail in the preceding 
sections of this chapter. This section will evaluate the 
procedures used to make the analysis. 

me Histribution Audit 

The distribution audit was a necessary first step 
in the data analysis. Although there was a multitude of 
information available concerning the cargo arriving at MOTBA, 
Only the number of shipments, their volume and their weight 
was of concern to this work. The audit took the raw data 
and placed it in a meaningful format. 

MircmailetribDUutton audit helped  tovestablish the fact 
that the Vietnam conflict generated a tremendous cargo flow. 
Mrerefore, it was decided not to use the six months of data 
ieeom FY ‘73. This left 52 weeks of data for FY ‘74, pro- 
frame a reasonable number of observations. Another alter- 
native would have been to exclude the data for the first 
two months of FY '74, because this data had large volume 
Mmaneriatilons on a scale similar to FY '73 data. The ex- 
elmision of this data may have resulted in the gencration of 
meer deviations for any or all of the forecasting techni- 
mies that were applied. 

Since there was a requirement to book container 
space aboard vessels three weeks in advance, the data was 
grouped in weekly intervals resulting in 52 observations 


mothe period. Although biweekly data could have been uscd, 
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it was more convenient to use weekly intervals for the data 
analysis. Intervals greater than two weeks would have re- 
sulted in too few observations. 

The distribution audit also identified the number 
of shipments classified as release unit because of weight. 
These shipments were considered part of the random volume 
flow based on the fact that they comprised a very small 
Bereentage of the total number of shipments. 

In short, the distribution audit provided a start- 
ing point for the procedures that followed. It simplified 
the data by placing it in meaningful perspective. 

fee oectial Correlation Test 

iivemoterewasealialyzed €@ reveal the presence of 
weet correlation. The Durbin-Watson d-Statistic Test is 
Sommoniy used for the purpose. The calculations for the 
test were easily accomplished by computer programming. The 
data was tested using weekly and biweekly intervals. Tri- 
Meckly intervals were not used since they would have re- 
sulted in only 17 observations. The results of the test 
feistrtied the application of the forecasting methods pre- 
viously discussed, 

Syme rorecasting Methods 
a. Moving Average Method 
The moving average was the most obvious method 
PoOminitially apply to the data. It provided a simple means 
of forecasting future volume inputs. Less weight was given 


to the weekly input as the time interval for the moving 


16 





average was increased. Time intervals from one to ten weeks 
were evaluated. To carry the evaluation beyond ten weeks 
would have been superfluous since time intervals beyond ten 
weeks result in minimal changes in the forecasts for each 
new input of data. As the moving average intervals become 
large, the process approaches the use of the mean of the 
mowume inputs over a certain time period (for example, a 
mmearter) as a forecast for the next similar period. The 
calculations indicated that increased smoothing of the data 
meseaccomplished by increasing the time interval of the 
moving average beyond three weeks. 

b. Exponentially Weighted Average Methods 

While previous discussions may have indicated 
Ghat there were four forecasting techniques applied to the 
Gata, in reality there were only two -- the moving average 
and the EWA plus two variations. 

The EWA method, unlike the moving average, con- 
Peers all past data to some extent. Through the applica- 
tion of the smoothing constant, the emphasis placed on the 
most recent data can be varied. 

One variation of the EWA was to adjust for trend 
effects in the data. The apparent trend was then smoothed 
and combined with a smoothed forecasted average to make a 
forecast for the following week. 

iiemsSececondavariationm was the application o£ an 
maaptive response rate to the EWA. This particular ARR used 


mesmoothed error term to calculate a tracking signal. Other 
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mecnods for calculating a tracking signal may be available; 
but they were not investigated. 

When the EWA methods were applied to the data, 
deviations were calculated beginning with the third week. 
For the moving average method, deviations were not computed 
Mem the week following the length of the time interval. 

In order tg evaluate the methods it may have been more ap- 
Meopriate to begin calculating the deviations for the EWA 
Meenods in the eighth, ninth or tenth week. 

4. Data 

The time period covered by the data was another 
memeeati1on in the application of the forecasting techniques. 
imeeeeneral, time-series analysis covers a three to four year 
period with data being grouped by months or even quarters. 
Conditions for this analysis were far from being optimal. 

The fact that there was only a year's data also pre- 
GCluded investigation of seasonal or cyclic effects. Again, 
meenree to four year period is usually required before 
smoothing techniques can be applied to the data. 

In summary, the breakdown of the yearly data into 
weekly intervals constituted the most logical approach to 
the problem. While a year of data did not represent a long 
@eriod of time, it did represent the best that was available 


ae the time. 
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V. CONCLUSION 


this thesis concerned itself with the forecast of the 
ferume Of cargo flowing into the Container Stuffing Station 
meeriolBA destined for a particular port of debarkation. 
iitesdata used in this work was taken from actual operations 
at the CSS during FY '74. Although every system is capable 
of generating errors, it was assumed that any errors were 
minimal and that the data was correct. 

Pater tO utilizing any forecasting techniques, the data 
meomanalyzed in three phases: a distribution audit; the 
BProuping of data into weekly intervals by POD; and a test 
Pemeserial correlation. 

iis thesis made use of the distribution audit of 
McCarthy and Carter. The purpose of the audit was to sim- 
Meaty the large volume of raw data generated by the opera- 
fens at MOTBA. The magnitude of the operations at the 
Seont’s illustrated by the following information. During 
fee /4 there were a total of 124,592 shipments weighing 
ever 176 million pounds and with a total volume of 10,070,724 
miipge feet. These shipments were destined for 3,356 consig- 
Mees in 157 overseas ports. Sixteen major ports accounted 
for 93% of the shipments, 93% of the volume and 91% of the 
weight. These 16 PODs then became the focal point for 
further analysis. 

iijemenouping of the data into weekly intervals by POD 


indicated the large volume fluctuations that occurred for 
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Meemeot the 16 PODs under consideration. Plots of the 
weekly volume inputs revealed the lack of any set pattern 
of cargo flow from POD to POD. Before making any decisions 
as to which forecasting techniques to apply, it was neces- 
sary to determine whether the data was random in nature or 
correlated from one observation to the next; that is, 
Serially correlated. The Durbin-Watson d-Statistic Test 
Satrorted this requirement. A computer program was written 
Eemcompute the d-statistic and to compare it with critical 
memees. Using weekly volume inputs the test indicated that 
fot !6 PODs had serial correlation in the data. Another 
Peat was performed by grouping the data into biweekly in- 
tervals. Similar results were obtained. The preceding 
results were then tested to determine if there was seriai 
Semnelation present in the data. The test was affirmative. 

The results of the serial correlation test supported 
the use of forecasting methods other than ordinary least- 
mares trending. The simplest and most logical forecast- 
meemethods to apply to the data were the moving average 
method and the exponentially weighted average method. Two 
Variations were applied to the exponentially weighted aver- 
meenresulting in a total of four methods. 

In order to evaluate the forecasting methods on an 
equal basis, the mean absolute deviation and the mean de- 
Viation squared were Geer EnCudeyvilattonepeine the ditfer- 
em@ee between the forecast and the actual input. Four 
computer programs were written to calculate the MADs and 


Mbeseter the various forecasting methods. 
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Moving averages from one to ten weeks were applied to 
the data. The MAD and MDS for each time interval were 
calculated by computer program. The EWA method used a 
smoothing constant to generate the forecast. Various levels 
femine Constant were applied to the data. A variation of 
the EWA calculated the apparent trend of the data. This 
apparent trend and the forecasted average were then smoothed 
and added to get the forecast for the following week. An- 
other variation to the EWA method was the application of an 
adaptive response rate. The ARR uses a tracking signal 
which is computed by taking the sum of the differences be- 
feen the forecast and the actual input volume for each 
week and dividing this total into the sum of the absolute 
Memes of those differences. As each error becomes avail- 
apie it is smoothed along with the old error. The smooth- 
mmemeonstant to be used in forecasting is then set equal to 
the absolute value of the tracking signal and then the EWA 
meeapplied as before. 

The deviation statistics resulting from the application 
of the above forecasting methods to the actual data with a 
one-week time horizon were tested using the Student t-dis- 
tribution. This test revealed that the EWA and ARR methods 
using the individual POD's optimal smoothing constants were 
Statistically the same and better than all other methods. 
Mires nbWwA technique was then chosen as preferred based on 
its computational ease relative to ARR. 

The selection of the EWA method offered some advantages 


over the three-week moving average being used at the CFD. 
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meest, EWA was computationally no more difficult. It de- 
rived its advantage by responding te changes in volume more 
uniformly than the abrupt alterations that the three-week 
moving average caused in the forecasts. Also, the accumula- 
tion of a historical data base in the smoothed forecast 
average increased the accuracy of the EWA method as the time 
horizon ere adsed. NeChe acc sarge 429 crm | Uo a tiGernO Ole. It 
mean absolute deviations for the one-, two-, and three-week 
time horizons using the EWA method evidenced the increase 
mamaecuracy. For this application and data, there can be 
little doubt that the use of the EWA provided improvement 
fmreche Lorecasting techniques. 

The other forecasting methods retained their relative 
merits. Greater accuracy was achieved using moving averages 
leeeencreasing the lengths of the averaging interval at the 
expense of increasing the computational difficulty. The 
apsence of any apparent trend in the data resulted in lack 
of any additional benefits from trend adjustments. While 
the adaptive response rate technique showed no real advan- 
tage over the EWA method for this data, its ability to 
respond rapidly to major policy changes makes it a viable 
Mieternative at the expense of increased computational dif- 
meeulty. As more data becomes available, all of these 
methods deserve additional consideration and thoughtful 
metention as the results of this work are based on only a 
Bear's data. 

eerowolh Sionificant increases in forecasting have been 


achieved, there remains in the data large amounts of randomness 
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Samsing significant errors. Its random nature, without 
jmeemas Or significant shifts, duave ued to the evolution of 
meme xponentially weighted average method as the preferred 
eaeece. As more data is accumulated the other methods may 
be reexamined. The optimal smoothing constants should be 
updated periodically to account for changes in the data and 
Meomote long run efficiency of the method. However, the 
most important safeguard against costly mistakes due to 
random effects in the input volume is careful managerial 
scrutiny of each forecast generated and used to make a 
booking. Judgment on the reasonableness of a prediction 
makes the manager the final and ultimate smoothing factor 


Mmimeene Utilization of the forecast. 
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APPENDIX A 


MEEK SHePMENTS VOLUME WEIGHT 
aa 3), 2UI8 220,992 By eon kO 
28 5 0d: AS 50 Or pe LOS 4S 
29 4,598 5 SOR a OO 7s 7 
30 4,884 SO S550 On 055 yrss 
BL Dey 344,109 O52 4502 
32 Amie eS 6,281,837 
oe 4,726 358, 880 OA ONO, FS 
34 ees 505. 9720 Socio 0 
55 Soe 378,598 1 Ol eps Os 
36 Seow S240 br 7 OG Od 
i) Bas 540447 6, 0041S 
38 4,638 501,952 5,808,959 
39 4,384 503755:5.0 Bo eZ 
40 Ae ya 300,089 Sle oa 
41 4,646 5555002 On 4 5 STAG 
42 aye oS Soro 0 O64 sey 
43 asi oS SOO oo 7 
44 Sy Se Comte aed Se 0 2S 
45 3,906 Moo o ko 4,840,206 
46 Se orks 248,060 AA Soo Ss 
47 570.0 0 288,476 5,085.827 
48 Sy 2 Sula S A205. 299 
49 oS 5204 048 5902, O07 
50 Seo D See ISS SIS 
ot S075 Sor 5 a5 Sec gooS 
52 Bao 79 BO ao Oe 2452 015.0 
Note: (1) Volume is in cubic feet 

iZeeenecight is ih pounds 
iaphe tl. Gross Weekly Inputs. Weeks 27 through 52, FY '73. 
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WEEK SHIPMENTS VOLUME WEIGHT 

1 2,444 261,448 S92 249 
2 25797 287,444 eee 0S 
5 eee 278,991 A ooo 54 
4 A WSs 252,296 Spear 42 
5 3,414 284,814 Sea Soo 
6 4,044 308,446 Do 2a 
iz 965 88,024 es 204 
8 ZOOS 208,179 Sy lorena 4 
9 Zoeell 196 ,694 SoZ omS 
10 1,994 oa 0 SUE OLS 
el 2,425 206% S55 SO aon © 
az Co ont DSc oer 
1.5 2,485 199,702 A oye 0) 
14 Ze W92, 8635 5,000 R 0g 7 
#5 1,784 So S06 Zo OO 
16 2655 Ia 7 Se one tS 
a 1,983 ones Ooo ole 
18 oo? i eo 2,199,669 
19 2,340 208,560 S766") 
20 Call 84,520 5, 2U0mb4 3 
zal Poo OS. O25 Oa oll) 
ZZ Deen 1495722 Sao on oe 
23 Zee 203,497 Soe uo 
24 Dod 199,888 S055 4025 
25 Zero 189 ,969 5 Ome 57 
26 LES 275 LES Ovo sles 2522 305 
Memes {1} Volume is in cubic feet 


(2) Weight is in pounds 


intitestl, Gross Weekly Inputs. Weeks 1 through 26, FY ‘74. 
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WEEK 


Zi 
28 
Zo 
30 
31 
32 
33 
34 
55 
36 
S/ 
58 
SS 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
a 
52 


Mette: (1) 
(2) 


ftaple II. 


See MeN ho 


1,530 
2,956 
3,058 
2,338 
2,580 
2,814 
2,667 
1,585 
3,007 
2,083 
1,994 
2,198 
2.557 
ee 
2,060 
2,145 
2,041 
2,704 
2,381 
2,189 
2,739 
1,395 
Oe 
Pesee 
2,812 
D2 8 


Weight is in pounds 


Gross Weekly Inputs. 
FContinued ) 
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VOLUME 


135,843 
232,446 
260,123 
179,400 
241,953 
242,661 
211,740 
139,665 
267 On 
164,209 
143, 360 
145,621 
172,236 
170,630 
201,845 
195,070 
POO 
198,084 
204,956 
193,582 
yates 
96,791 
185,867 
194,153 
229,928 
188,520 


Volume 1S “In Cubic feet 


WEIGHT 


2,186,285 
ease 
4,468,724 
3,078,738 
Beis 19 
4,034,761 
EO aes 
2,459,631 
4,907,835 
D725), 845 
2,385,097 
DE ia 
3,151,837 
2,921, 166 
3,388,849 
3,211,454 
3,135,205 
3,052.520 
3,269,329 
3,264,435 
SSIS 7.0 
1,628,575 
ee esol 
2,716,456 
3,602,602 
Poole 60s 


Neco? stirough) S25 FY 


PAs, 





Period analyzed - FY '74 365 days 


Shipments 7459 2a sh pmenes 
Volume of cargo ORO? 07 24 el wes reet 
Weight of cargo Or 1 5.555 pounds 
Number of PODs oe 

Number of consignees Sooo 


maple Tit. Data Summary for the Distribution Audit. 
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SHIPMENTS 


Total shipments of 16 PODs LAS ees 
Total shipments of all PODs 124,592 
Percentage of total shipments 92.89% 
VOLUME 
Total volume of 16 PODs Soe On Ors 
Total volume of all PODs TORO Ge 24 
Average volume/shipments (16 PODs) 87 
memeentage of total volume (16 PODs) OS 055 
WEIGHT 
Total weight of 16 PODs HOO SS 7 S272.6 
Total weight of all PODs LO ea Se SoS 
Average weight/shipments (16 PODs) eo ore 
Percentage of total weights (16 PODs) 91% 


shipments 


shipments 


cubic feet 
cubic feet 


cubic feet 


pounds 
pounds 


47 pounds 


igbae iy, Tabular Data for 16 Major PODs, FY ‘74. 
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POD 


RGU 
RG1 
RJ1 
RJ 3 
SA3 
TAI 
TA2 
UB1 
UC2 
UD6 
UL7 
UM1 


XE2 
XJ1 


Note: 


SHIPMENTS VOLUME 
20 AOE 1,474,747 
6,041 491,465 
S50 iO OIG 
2,848 170,635 
Za2 LO 106,150 
16,837 1,377,838 
2,502 2 OmOn 
6,793 502,989 
12,339 1,081,519 
6,593 976,869 
8,691 451,466 
3,183 msn 
6,619 434,615 
3,769 276,100 
14,008 1 Sls 559 
1,324 186 ,690 


Cimevolume iS in Cubic treet 


(2) Weight is in pounds 


Taple. Vv . 


oor 


POD a tae Oe LO Vice 


WEIGHT 
28,284,704 
11,486,615 

3,774,348 
BS 255555 
eeareOt2 
19,920,384 
Delisi2 e227 
7,269,693 
19,311, 360 
19,325,856 
9,062,837 
4,445,639 
7,081,348 
3,740,452 
16,679,435 


2,935,428 





NUMBER OF PERCENTAGE SOF 


POD RU SHIPMENTS RU SHIPMENTS 
RA3 384 1.90 
RGU 156 2.58 
RG1 37 2.00 
RJ1 41 es 
RJ3 15 0.68 
SA3 183 1.00 
TAL 15 0.59 
TA2 86 1.26 
UB] 242 1.96 
UC2 B27 4.90 
UD6 80 0.92 
UL7 43 1.35 
UM1 80 0.92 
UM4 20 0.53 
XE2 136 0.97 
ed 55 4.15 
TOTAL 1,900 


lime percentage of RU shipments to the’ total number of 
shipments for the 16 major PODs is 1.64%. 


Table VI. Release Unit Shipments. 
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Week WEEK NO. d-STATISTIC BEEN WHEKe® NO@rd-STATISTIC 

Eee INT START OBS MEANING EOm Nl START OBS MEANING 
RAZ 1 i a INCN UB1 1 1 a2 NONE 

2 1 26 INCN 2 if 26 NONE 

2 2 aS POS 2 2 DS NONE 
RGU 1 1 52 POee UGZ.-~ vl 1 Be NONE 

Z af aS INCN 2 I 26 INCN 

2 2 ZS INCN 2 2 ZS INCN 
RG1 1 i 52 Poa UD6 1 1 52 INCN 

2 1 26 POS 2 1 26 NONE 

Z 2 25 Peo 2 2 25 NONE 
Roe 1 1 52 E@ol UL7 1 i BZ POs 

2 il 26 NONE Z if 26 Foe 

2 2 25 NONE Z 2 25 POG 
RJ3 1 i} 52 NONE UM1 1 1 52 Poot 

2 il 26 NONE Z if 26 NONE 

2 Z ZS NONE 2 Z 25 POs 
SAS 1 if oe POSI UM4 1 1 52 POSI 

2, 1 26 INCN Z i! 26 NONE 

2 2 as POSI 2 2 25 EOow 
eal 61 i! ne INCN Xs Zee 1 52 POST 

2 1 26 POs 2 if! 26 POS 

2 2 Zoe INCN 2 2 25 Rost 
a2 41 i 5Z INCN ee 1 52 EOow 

2 i 26 INCN 2 il 26 Posi 

2 2 25 INCN 2 Z ZS INCN 
femee, (1) POSI = the test indicated POSITIVE serial correlation 

(ZI NGN = the test was INCONCLUSIVE 
ComecNre= the test @idgcaced no Serial correlation 
iaes vil, Serial Correlation Test Results. 
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Table 


POD TWMA 
RA3 6797* 
RGU 3910 
RG1 1754 
RJ1 1793 
RJ 3 1130 
SA3 7862 
TAL 996 
TA2 3654 
UB1 B52 1 
UC2 5 56 
UD6 4205 
UL7 1383 
UM1 Boo, 
UM4 Zoe I 
MEZ 7770 
me 1 ae 5 
Note: (1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(eet) 
*Indicates 


{ct 


> 


10 


10 


10 


MAD 
POD 


TWMA 


OMA 
t = 
EWA 
Qa = 


EWA- 


EWA- 


y= 


VIIT3? MADS} forme Pels. 
One-week time horizon 
OMA EWA ‘EWA-TA EWA-ARR 
3 2 ate 

Oo ee Ome) Cole eco = 78. 7074 
Oy) One Leo Siac 5 MSL 7 240 5050 
EAC eee Oley OO 547 87 soa” *1GU1 
G2 es Oe OA ey 7. OMe lei ard 
od 7 eS 5, eno BOI 1G s0- 
(esr ee ey hy eee) SOnlen eign Oh7 4 

ol tee22. 955 @7..0) Jol / 1S aos 
OS ted Cleon oe OunoaOI = iw 4.04 
SOO tl et Omori 4 = 55 5 216 
a5c5” 2 lhOrao40 W044 748g9.- 17 4794 
Soe COO Oa meoro a7 ae S57 oo 
ZG, =. 210 i290 er, 2 alee 8 OOM ecw? 
Od es Se 0 ee Oars neler 0 eS Lol 
7500 ee 40 roo Ce eee lee OSo% 
15947 27 6 See Oey now Oe eee 510 
2510 1 Ses oe Oran OO) 2475 
= Mean Absolute Deviation 
=) Orte on Nebarkacronm 

= Three Week Moving Average 
= Optimal Moving Average 
time interval (an weeks) for OMA 
= Exponentially Weighted Average 
SNOOCHIN co COns Galea cumENurand EWA-TA 
TA = Exponentially Weighted Average with 

Trend Adjustment 

ARR = Exponentially Weighted Average with 


Adaptive Response Rate 
smoothing constant for EWA-ARR 


MADS are expressed Pn cuble Leet 


smallest MAD for POD 


7 





apie VIII. 


POD WMA 
RA3 6204 
RGU SoZ | 
RG1 1601 
Ro i Nesi5 0 
Rod 3 850 
SA3 6639 
TAI] 798 
TA2 3194 
UBL 4655 
Wig 2 4001 
UD6 5688 
Wie 7 1053 
UM1 3146 
UM4 2343 
XE2 6984 
Pe) 1 we] |. 


MADs for 16 PODs (Continued). 


Two-week time horizon 


+ OMA EWA q EWA-TA | EWA-ARR 
3 6204 .32 5955* .17 6293 - .89 6095 
3 3621 .09 32.84*.22 3489 .45 3339 
tO eles eS nnS)7 GOO .23 1311* 
8 1194* .37 1420 .17 1449 50) 128% 
S50 (a3 80 bl Aid .01 750% 
S 6659 402 6250 UL G2 5 VS aoey 
6 95e5* 622 Fone .12 729 .11 694 
10 3056 .39 3072 .17 3125 - .03 2991% 
7 4250 .15 4246 .06 4179* - .01 4417 
7 3150* .15 3399 .06 3278 - .10 3423 
723155 427 3285 .l0eszo8 = 608 2938 
4 1005* .28 1013 .17 1082 - .44 1033 
10 2641 .19 2816 .09 2893 .01 2615 
9 2153 -.02 1878 -.01 1847* 0 ey 
36984 .65 6697 .31 6941 - .01 6269* 
1D D2 2S 2085. 10 2590 “> Sey, 


Emrdieates Smallest MAD for POD 
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Table VIII. 


POD 


RGU 
RG1 
RJ1 
RJ 3 
SA3 
TAI 
TA2 
UBL 
UC2 
UD6 
UL7 
UM1 
UM4 
pee 2 


XJ1 


eimdicates 


TWMA 


6620 
5477 
SZ 


1188 


6451 
ial 
3149 
io. 7 
5 SO 
Sr 210 
De 2 
O25 
21 
oro dl 


| 


MADs for 16 PODs (Continued). 


| c+ 


10 


10 


10 


0 
10 


10 


10 


EO 


Ematllest MAD ter LOD 


OMA 


OReo, 
3410 
0p 


i Sito 


6388 
545- 
Too s 
ore 
BS AS © 
Aiea 
oes 
Loo 
970 
ply 


2049% 


aS 


ZO 
li 
eS 
43 
Sone 
~ 34 
oO 
ork 
pn 
mao 
eo 
cS 
eg 
Oe 
OD 
ld 


EWA 


Sona 
Saez Oe 
1498 
EZ55 
Og 
5848 
626 
3030 
S205 
Zo a0 
Zs 
946 
ale) (5 
aoe 
G7 7S 


Zou s 
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Three-week time horizon 


EWA-TA 

0 
-10 6069 
pO emo Ore 
5 IES Ae 
Poe 1 296 
.09 784 
POL a 2a0. 
13 653 
ie 2990 
Oe OHO = 
OG we cyecn 
oy 
yGe OOS 
UO 2 oed 
Oe 20a 
Toe OUG0 
OST ae 


=e 


EWA-ARR 

Ne 
poy Oot 
47 3214 
coo) eae 5% 
One AO.8 © 
2 Oi halos 
00 6234 
14 S071 
B0S9 62759 
OUTS Sa a 
UOT ee 20-512 
rO7 ee cao 
46 Oe 
elec, 5 onl: 
eae Ao, 
Odeo 7 
Oc oO 





EOD 


RGU 
RG1 
RJ1 
RJ 3 
SA3 
TAI 
TA2 
UB1 
UG2 
UD6 
UL7 
UM1 
UM4 
pele 2 


XJ 1 


Note: 


TWMA OMA EWA EWA-TA EWA-ARR 
t = e at 
oa COUN Z 50. ee 60) 958540206 98,507 = 182. 985528 
mie > 6 427 546 oe 7 ee a 20090 woo (255,060 
10,349 10 5,769* .60 10,736 .30 10,117 .25 #&8,805 
wees Ge691* 335 965.337°1.16 5.511 .02 5,474 
Ueoce. f 255991 gdeZ PAS See Oop eed Oo ly Fag So ZS S 
Domo See Veet SO" 4 56am 96, /2ouencd 104,177 -1.00 104,207 
coco el isd 2) 545 io 611°. 11 9 1,614 
eo se 0 Sao Ss 28 52054565 12 215253 = .035- 20,549 
51,088 ANOS? Wiel Tee OO elon 4.5547 9-76 44, 860 
41,418 30,695* .14 35,775 .05 34,811 -1.00 35,382 
pe 996 SOPUMiCn Ney oe OS (emo ol 2017 07 ey 5 7 OU 
eons Sol OS Aas Se Omens 5.20.71 -~ 474 Sloe 
Sel Oso 0S 5 ee 2 ee Omer 5 0105 =" 145 21,4637 
elo 7 O OPA Oo aes Ts Omen o LUN = ee So 7105500 
Miers 10 84°5546* .$271035,643 9.25 106,797 OS Oo 5 2 
15,387 135 SO Semi ae ome 1 1S 5 2, 400 
(1) MDS = Mean Deviation Squared 
(2) POD = Port of Debarkation 
(3) TWMA = Three Week Moving Average 
(4) OMA = Optimal Moving Average 
(jee - time interval (insweeksj@ror OMA 
(6) EWA = Exponentially Weighted Average 
(7) a = smoothing constant for EWA and EWA-TA 
(8) EWA-TA = Exponentially Weighted Average with 
Trend Adjustment 
(9) EWA-ARR = Exponentially Weighted Average with 
oapeiye eke Semis ce wake 
(10) y = smoothing constant for EWA-ARR 
(11) MDSs are in 1,000's and expressed in (cubic feet)? 


Table IX. 


MDSs for 16 PODs 


One-week time horizon 


mumadweates ,smallest MDS for POD 
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POD 


RGU 
RG1 
RJ1 
RJ3 
SA3 
TAI 
TA2 
UB] 
UC2 
UD6 
UL7 
UM1 
UM4 
1 OA 


XJ 1 


TWMA 


70,264 
oh S34 
9,064 
B52 
1,659 
76,062 
1,060 
16,801 
30,872 
25,069 
24,416 
2,084 
19,208 
10,062 
81,399 


oe 7 3 


apace Lx. 


| ct 


10 


10 


10 


10 


10 


10 


oe) 


10 


MDSs for 16 PODs 


(continued) 


Two-week time horizon 


OMA 


69,626 
KK 
5 293% 
Oy Asieg 
1,388 
76,062 
749% 
Te 
24,900% 
14 ,927% 
20,683 
1,880% 
15,958 
6,781% 
80,557 


10 ,908* 


2 


ne 
meer 
69 
Se 
oS 
07 
bee 
226 
14 


ee) 


a0) 
ae 
als 
.o4 


so 


EWA 


62,116% 
18 ,672% 
8,385 
4,019 
1,327% 
73,856% 

953 
14,376 
25,001 
18,843 
Dianova 

noice 
16,063 
7,909 
79,063 


11,087 


Mimdicates smallest MDS for POD 


uo 


[R 


205 
snilsd 
150 
kG 
2 V6 
705 
we 
Alba 
UG 
. 06 
als 
cae 
aalt@) 


see 


ee 


EWA-TA 


64,218 
iste 
Teele I 
4,165 
1, 360 
73,992 
eons 
15,107 
25,749 
17,843 
22,598 
Ze Oa 
16,948 
8,057 
81,152 


ea 6 


Bore 
eS. 
nee 
0 
oD 
.00 
et 
UES 
aaa 
.00 
aes 
47 
eel 
mel 
ua 


14 


EWA-ARR 


69,925 
18,829 
5,588 
3,965 
1,334 
80,384 
948 
es aon 
Ber ase 
17,496 
17 ,668% 
eee 
14,948% 
8,216 
60,822% 


11,067 





POD 


RGU 
RG1 
RJ1 
RJ 3 
SA3 
TA1 
TAZ 
UB1 
UC2 
UD6 
UL7 
UM1 
UM4 
OE Z 
XJ1 


TWMA 


On OZ 


19 ,006 


10,187 
PTS 7 
1,347 

74,786 

806 

15,444 

18,965 

OS 

20,080 
1,633 

17,862 
787 

81,188 


isn. 778 


Table ie 


Gat: 


10 


10 


10 


10 


10 


HO 


10 


10 


10 


MDSs for 16 PODs 


(continued) 


Three-week time horizon 


OMA 


64,618 
18,391 
5, 461% 
pele 
1,035 
74,786 
663 
9,783% 
16,064 
10,170% 
15,516 
1,551% 
12,662 
5,197% 
75,150 


SA 106% 


*Indicates smallest MDS 


[2 


Si) 
220 
ais 
noo 
14 
US 
a6 
24 
ee 
. 16 
24 
OO 
io 
sun 
oO 
a 


ions 


ov 


EWA 


54,579% 
16 ,418% 
8,571 
3,234 

967% 
61,446 
7A 
eS 

15,866% 
13,467 
16,499 
1,570 
esi 
6,067 
79,907 


10,803 


POD 


[2 


04 
2 ek 
oe 
ks 
76 
WZ 
14 
a0 
Oy 
.06 
rio 
mes 
pale) 
.O1 
20 
oz 


EWA-TA 


55,583 
18,245 
6,724 
3,354 
1,002 
61, 394% 

779 
13,408 
16,757 
12,539 
18,160 

1,638 
14,584 
6,217 
80,197 


le 20 


78 
ONL 
24 
rl 
me) 


lt 
NS 
mut 
.00 
U7 
249 
BZ 
HOG 
peel 
94 


EWA-ARR 


67,596 
16,570 
5,822 
2,968 
969 
73,494 
635% 
10,545 
17,603 
11,723 
12,549% 
L211 
11,984% 
7,084 
54,919% 


10,820 
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pue YM JOF JYUueYSUOD BuTYAoOOUS = VD (6) 
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WEEKLY 
TIME TWMA OMA EWA EWA-TA ARR 


HORIZON ic O O Y 
one — Soa ome. 1 Oo See One 57? OSS, 979 .17 566,904 -.14 $43,745 
two MO Jo Omen ue 469 -156 7.560 ,940 417 421,452 =.05 393,060 


three 375,569 10 341,400 .36 350,451 .17 384,810 -.05 346,852 


Note: 
(1) MDS = Mean Deviation Squared 
(2) OMA = Optimal Moving Average (a specific optimal interval 


Mmemapplied to all POD data) 

(3) EWA = Exponentially Weighted Average (a specific optimal 
fee vel is applied to all POD data 

(4) EWA-TA = Exponentially Weighted Average with Trend Ad- 
justment (a specific a-level is applied to all POD data) 

me) ARR = Adaptive Response Rate (a specific optimal y-level 
iomapplied to all data 

(6 ) 


t = optimal time interval for OMA 
(7) a = smoothing constant for EWA and EWA-TA 
(8) y = smoothing constant for ARR 


(9) MDSs are expressed in (cubic feet)? 


Table XI. Aggregate MDSs for 16 PODs. 
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imethis work the Student t-distribution was used to 
resolve whether the mean of the MAD differences between two 
forecasting methods was equal to zero. 

mie hypotheses were: 


H,: l= Ou. H,: u # OQ. 


mae t-test utilizes the following 


oe (X - HM) Yn 


where X is the sample mean, w=0, n=number of observations 


and s=the sample standard deviation. 


Pabculated values of t are compared with critical values 


Cit. If 
ta(eauke Ulbaeech) << tose rd timed | ) 
H. is accepted 
If t (ealeulated), > tert mueal) 


Ho PSeNOt accepecd: 


Tt Le Poet acceptor Evommeumods are qdurterent and 
the method that generates the lowest mean absolute deviation 


[emenosen as the preferred. 


Mitemtest only Wtilized deviations generated by a one- 


week time horizon. 


fabless ot -lesturrocedure. 


100 





COMPUTER OUTPUT ONE 


TABLE OF STAFISTICAL DATA FOR SFLECTCO POUS USING WEEKLY TOTAL VOLUME (CU. 


FT.) DATA 


OHS. MEANING 
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COMPUTER OUTPUT THREE-ONE 


Mean Absolute Deviation (in cubic feet) 
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COMPUTER OUTPUT THREE-THREE 
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COMPUTER OUTPUT FOUR-ONE 


Mean Absolute Deviation (in cubic feet) 
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COMPUTER OUTPUT FOUR-THREE 
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COMPUTER OUTPUT FIVE-ONE 


Mean Absolute Deviation (in cubic feet) 
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COMPUTER OUTPUT FIVE-THREE 
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COMPUTER OUTPUT FIVE-FOUR 
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COMPUTER PROGRAM ONE 
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TICS AND MEANING FOR WEEKLY INTERVALS. 
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-STAT FOR SAMPLE WITH 52 OBSERVATIONS. 


RMINING MEANING OF D-STAT FOR SAMPLE WITH 26 OBSERVATIONS. 
RMINING MEANING OF D-STAT FOR SAMPLE WITH 25 OBSERVATIONS. 
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